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ABSTRACT  We  present  experimental  and  theoretical  results  on 
white-light  generation  in  the  filamentation  of  a  high-power  fem¬ 
tosecond  laser  pulse  in  water  and  atmospheric  air.  We  have 
shown  that  the  high  spatio-temporal  localization  of  the  light 
field  in  the  filament,  which  enables  the  supercontinuum  gen¬ 
eration,  is  sustained  due  to  the  dynamic  transformation  of  the 
light  field  on  the  whole  transverse  scale  of  the  beam,  includ¬ 
ing  its  edges.  We  found  that  the  sources  of  the  supercontinuum 
blue  wing  are  in  the  rings,  surrounding  the  filament,  as  well 
as  at  the  back  of  the  pulse,  where  shock-wave  formation  en¬ 
hanced  by  self- steepening  takes  place.  We  report  on  the  first 
observation  and  demonstration  of  the  interference  of  the  super¬ 
continuum  spectral  components  arising  in  the  course  of  multiple 
filamentation  in  a  terawatt  laser  pulse.  We  demonstrate  that  the 
conversion  efficiency  of  an  initially  narrow  laser  pulse  spectrum 
into  the  supercontinuum  depends  on  the  length  of  the  filament 
with  high  intensity  gradients  and  can  be  increased  by  introduc¬ 
ing  an  initial  chirp. 

PACS  42.65.Jx;  42.65.Re;  42.25 .Bs;  42.50.Hz 


1  Introduction 

Propagation  of  a  high-peak-power  ultra-short  laser 
pulse  through  bulk  transparent  media  is  accompanied  by  the 
strong  modification  of  its  spatial  and  temporal  properties. 
Transformation  of  the  pulse  into  the  frequency  domain,  in  par¬ 
ticular  the  appearance  of  a  large  frequency  sweep,  which  can 
extend  from  the  ultraviolet  to  the  infrared  range  (supercon¬ 
tinuum  generation  [1]  or  a  white-light  laser  pulse  [2]),  is  the 
spectral  representation  of  the  spatio-temporal  modification  of 
the  radiation  in  the  nonlinear  medium. 

Earlier  experiments  on  the  nonlinear  pulse  transformation 
and  broadening  of  the  frequency  spectrum  were  performed 
with  nano-  and  picosecond  pulses  focused  in  condensed  me¬ 
dia  [3].  Advancement  of  laser  technology  creating  femtosec¬ 
ond  pulses  with  peak  power  up  to  many  terawatts  made  pos¬ 
sible  the  observation  of  the  nonlinear  pulse  shape  and  spectral 


Kl  Fax:  +7-095/939-3113,  E-mail:  kosareva2003@mail.ru 


transformation  in  gaseous  media.  One  of  the  first  observations 
of  spectral  superbroadening  from  ~  400  to  ~  800  nm  in  gases 
was  performed  in  [4],  where  70-fs  and  2-ps  dye  laser  pulses 
centered  at  600  nm  were  focused  in  Xe  and  N2  gases  at  30-atm 
pressure. 

Self-transformation  of  the  pulse  shape  and  spectral  broad¬ 
ening  are  the  result  of  strong  nonlinear-optical  interaction 
of  the  light  field  with  the  medium,  which  takes  place  in  the 
conditions  of  high  localization  of  the  radiation  both  in  space 
and  time.  There  might  be  several  ways  to  localize  the  light 
field  inside  the  medium.  One  possible  way  is  to  take  advan¬ 
tage  of  an  optical  device.  With  such  devices  as  capillaries, 
microstructured  fibers,  and  photonic  crystals  it  is  possible  to 
govern  the  mode  dispersion  and  to  achieve  high  localization  of 
the  light  field  at  comparatively  long  interaction  distances  (see 
e.g.  [5-7]).  Akimov  et  al.  [8]  reported  the  first  observation  of 
the  spectral  superbroadening  (up  to  400  nm)  of  an  800-nm, 
40-fs  pulse  with  subnanojoule  energy  in  a  tapered  fiber.  Thus, 
one  of  the  important  advantages  of  the  supercontinuum  gen¬ 
eration  in  microstructured  fibers  and  photonic  crystals  is  the 
relatively  low  input  pulse  energy  required  for  the  registration 
of  the  phenomenon. 

In  bulk  media  the  localization  of  the  light  field  is  achieved 
due  to  geometrical  focusing,  self-focusing,  or  temporal  com¬ 
pression  of  the  radiation.  This  localization  may  persist  for 
long  distances.  In  the  first  experiments  on  the  propagation 
in  air  of  pulses  generated  by  Ti:  sapphire  laser  amplification 
systems  with  peak  powers  of  5-50  GW,  duration  100-250  fs, 
central  wavelength  700-800  nm,  and  collimated  input  beam 
intensity  distribution,  the  phenomenon  of  filamentation  was 
observed.  In  the  course  of  filamentation  9%-12%  of  the  input 
pulse  energy  was  localized  in  the  filament  -  a  narrow  near- axis 
region  with  the  transverse  diameter  of  the  order  of  100  \im  and 
a  longitudinal  scale  of  several  tens  of  meters  [9-11].  Later 
on,  a  more  than  two-hundred-meter  filament  was  obtained 
in  the  propagation  of  a  60-fs  pulse  centered  at  795  nm  with 
a  peak  power  of  300  GW  [12].  With  a  terawatt  peak  power  of 
the  Ti: sapphire  laser  system  and  a  120-fs  pulse  duration  the 
filaments  and  accompanying  spectral  broadening  were  regis¬ 
tered  not  only  in  air  but  also  in  the  atmospheric  density  rare 
gases  [13].  There  were  observations  of  filamentation  at  the 
laser  central  wavelength  different  from  790-800  nm.  In  par- 
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ticular,  in  [12]  the  experiment  was  also  performed  with  1-ps 
pulses  at  1053  nm  and  in  [14]  450-fs  pulses  at  248  nm  were 
used.  The  analysis  shows  that  the  extended  localization  of  en¬ 
ergy,  i.e.  filamentation,  is  the  result  of  dynamic  competition 
between  the  Kerr  self-focusing  of  the  radiation  and  its  nonlin¬ 
ear  defocusing  in  the  self-produced  laser  plasma  [1 1, 15-17]. 
The  length  of  the  nonlinear  interaction  in  the  filament  can  be 
long  enough  to  transform  the  essential  part  of  the  input  pulse 
energy  into  a  wide  range  of  spectral  frequencies.  An  important 
advantage  of  supercontinuum  generation  in  bulk  media  is  the 
possibility  to  obtain  a  wide  spectrum  inside  the  medium  under 
investigation  (e.g.  in  atmospheric  air  or  oceanic  water)  and  to 
deliver  the  resulting  coherent  white  light,  or  a  white-light  laser 
pulse  [2],  to  the  prescribed  position,  which  can  be  at  a  very 
large  distance  from  the  laser  source. 

The  frequency  spectrum  produced  in  the  course  of  a  long- 
range  filamentation  of  high-power  femtosecond  laser  pulses 
in  the  atmosphere  incorporates  nearly  all  typical  features  of 
the  supercontinuum,  which  will  be  discussed  in  more  detail  in 
Sect.  2.  These  features  are  a  large  blue  wing  produced  mainly 
due  to  free-electron  generation,  a  red  wing  created  due  to 
the  combined  effect  of  instantaneous  and  delayed  responses 
of  the  Kerr  nonlinearity,  and  an  angular  distribution  of  spec¬ 
tral  wavelengths  in  the  form  of  conical  emission.  A  spectrum 
that  extends  from  300  nm  to  4.5  |im  after  the  propagation  of 
both  35-fs,  60-mJ  and  100-fs,  100-mJ  laser  pulses  centered  at 
800  nm  was  obtained  in  [18]  in  atmospheric  air.  In  the  first 
femtosecond  lidar  experiments  [19, 20],  time-resolved  broad¬ 
band  spectroscopy  was  performed  on  the  supercontinuum 
backscattered  from  the  atmosphere  from  the  height  of  more 
than  1000  m.  High-resolution  spectra  of  the  range-resolved  at¬ 
mospheric  transmission  showed  the  ro- vibrational  H2O  and 
molecular  oxygen  absorption  transitions.  Thus,  a  coherent 
white  light,  accompanying  filamentation  in  air,  proved  to  be 
a  novel  light  source  for  atmospheric  remote  sensing. 

The  purpose  of  this  paper  is  to  present  our  experimental 
and  theoretical  studies  of  the  self-transformation  of  a  fem¬ 
tosecond  laser  pulse  into  a  white-light  laser  pulse  in  bulk 
condensed  and  gaseous  media  by  the  example  of  water  and 
atmospheric  air,  respectively.  We  consider  the  joint  effect 
of  such  physical  mechanisms  as  Kerr  nonlinearity,  ioniza¬ 
tion,  group-velocity  dispersion  as  well  as  initial  pulse  pa¬ 
rameters  such  as  wavefront  curvature  and  frequency  chirp 
on  the  nonlinear  spatio-temporal  transformation  and  the  cor¬ 
responding  frequency-angular  spectrum  of  the  pulse.  By  the 
frequency-angular  spectrum  of  the  pulse  we  mean  spatio- 
temporal  Fourier  transform  of  the  complex  light  field  ampli¬ 
tude  (see  (15)  in  Sect.  4  of  this  paper). 

2  Mechanisms  for  supercontinuum  generation 

Let  us  consider  possible  mechanisms  for  the  for¬ 
mation  of  the  new  frequencies  in  the  continuum  spectrum  of 
a  self-transformed  laser  pulse.  By  the  term  ‘supercontinuum 
generation’  we  understand  the  continuum  range  of  all  the  new 
frequencies,  which  arise  due  to  the  nonlinear  transformation 
of  the  pulse  in  the  medium  and  we  exclude  the  spectrum  of 
the  initial  laser  pulse.  At  the  same  time,  the  term  ‘white-light 
laser  pulse’  encompasses  the  whole  spatio-temporal  light- 
held  distribution  and,  consequently,  the  modified  frequency 


spectrum  that  was  created  after  the  propagation  in  the  non¬ 
linear  medium.  This  total  light-held  distribution  also  includes 
the  spectrum  of  the  initial  laser  pulse. 

One  of  the  hrst  explanations  of  ultra-short  laser  pulse 
spectrum  transformation  was  done  on  the  basis  of  a  simple 
model,  where  the  maximum  relative  broadening  Aco/coo  was 
directly  proportional  to  the  pulse  intensity  in  the  medium  with 
instantaneous  Kerr  nonlinearity  [21].  The  frequency  deviation 
8co  (t)  =  —kzdAn/dr  (here  An  is  the  nonlinear  contribution 
to  the  refractive  index,  t  is  the  pulse  duration,  z  is  the  propa¬ 
gation  coordinate,  and  k  is  the  wave  number)  is  symmetric 
relative  to  the  peak  of  the  pulse.  Therefore,  the  spectral  broad¬ 
ening  is  also  symmetric  relative  to  the  laser  central  frequency. 
In  this  model  the  radiation  is  considered  to  be  a  plane  wave 
with  a  given  temporal  distribution.  The  spatial  distribution  of 
the  radiation  in  the  transverse  direction  was  not  taken  into  ac¬ 
count.  Later  on,  it  was  concluded  that  self-focusing  played  an 
essential  role  in  spectral  superbroadening  [4, 22, 23].  Indeed, 
the  simple  self-phase-modulation  model,  which  neglects  the 
spatial  effects,  could  not  explain  the  increase  of  the  spectral 
width  A to  up  to  several  hundred  nanometers,  observed  in  the 
experiments. 

Further  experimental  studies  have  shown  that  not  only 
self-focusing  but  also  geometrical  focusing  of  the  pulse 
into  the  medium  influence  the  width  of  the  frequency  spec¬ 
trum  [24].  It  was  demonstrated  that  the  threshold  of  the 
supercontinuum  produced  by  a  1-ps,  593-nm  pulse  focused  in 
20-atm  CO2  gas  increases  with  the  decrease  in  the  geometrical 
focal  length  of  the  lens.  This  threshold  increase  was  associ¬ 
ated  with  the  defocusing  of  the  radiation  by  the  plasma,  pro¬ 
duced  in  the  vicinity  of  the  geometrical  focus.  The  shorter  the 
focal  length  of  the  lens,  the  stronger  was  the  scattering  of  the 
radiation  by  the  plasma.  The  same  conclusion  was  obtained 
in  recent  experiments  using  45-fs,  800-nm  laser  pulses  on  the 
competition  between  filamentation  and  multiphoton/tunnel 
ionization  (in  air  [25])  and  optical  breakdown  (in  condensed 
matter,  water  [51,61],  and  fused  silica  [26]). 

At  the  same  time,  the  free-electron  generation  contributes 
essentially  to  the  spectral  transformation.  The  frequency  devi¬ 
ation  Sen  (t)  =  —kzdAn/dz  is  positive  due  to  the  negative  con¬ 
tribution  to  the  refractive  index  An  from  the  free  electrons,  the 
density  of  which  grows  during  the  pulse.  The  spectral  broad¬ 
ening  is  asymmetric  and  blue  shifted  relatively  to  the  laser 
central  frequency  [27-30].  Ionization  and  optical  breakdown 
(the  latter  often  follows  the  ionization  in  dense  gases  or  con¬ 
densed  media)  strongly  depend  on  the  pulse  peak  intensity. 

In  the  conditions  of  tight  focusing  of  the  radiation  into 
the  medium,  where  the  third-order  nonlinear  susceptibility 
is  comparatively  low  (e.g.  in  the  noble  gases),  the  major 
contribution  to  the  spectral  transformation  comes  from  the 
free  electrons.  Experimentally,  this  phenomenon  was  studied 
in  [28, 29]  and  numerically  in  [29-31].  The  extension  of  the 
blue  wing  formed  after  the  tight  focusing  of  a  100-fs  pulse 
with  the  central  wavelength  of  620  nm  into  the  noble  gases 
varied  from  ~  15  nm  in  1-atm  argon,  krypton,  and  xenon 
to  ~  50  nm  when  these  gases  were  at  5-atm  pressure.  Spec¬ 
tral  blue  shifting  is  accompanied  by  the  spatial  defocusing 
of  the  radiation  and  the  formation  of  rings  in  the  trailing 
part  of  the  pulse  [32].  It  was  found  as  the  result  of  numer¬ 
ical  simulations  [31,33]  and  registered  experimentally  [34, 
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35]  that  the  blue-shifted  components  of  the  frequency  spec¬ 
trum  are  formed  at  the  off-axis  position  in  the  transverse 
distribution  of  the  radiation.  A  larger  frequency  shift  is  gen¬ 
erated  at  positions  farther  from  the  propagation  axis.  Thus, 
the  spectrum  of  the  pulse  is  transformed  not  only  in  the  fre¬ 
quency  domain  but  also  in  the  domain  of  the  spatial  wave 
numbers. 

Spatial  characteristics  of  the  white  light  generated  in  con¬ 
densed  media,  in  particular  in  H2O  and  D2O,  were  studied 
in  [36].  The  supercontinuum  was  produced  by  focusing 
a  22-ps  pulse  at  1.06  |xm  of  up  to  35-mJ  energy  into  water. 
Independently  of  the  focal  length  of  the  lens  used  to  focus 
the  radiation  into  the  cell,  the  supercontinuum  was  spread  out 
in  a  circular  rainbow  with  frequency  increasing  with  the  in¬ 
crease  in  the  off-axis  angle.  Golub  [36]  interprets  the  observed 
phenomenon  -  the  supercontinuum  cone  emission  -  as  the  ra¬ 
diation  emitted  from  the  surface  of  self-trapped  filaments  at 
angles  that  satisfy  the  Cerenkov  relation. 

Another  possible  explanation  of  the  frequency  angular 
distribution  of  the  supercontinuum  was  suggested  in  [37], 
where  the  increase  of  the  off-axis  angle  with  increasing  spec¬ 
tral  frequency  of  a  100-fs,  l-|xJ  pulse  focused  into  ethylene 
glycol  was  attributed  to  four-photon  parametric  generation  on 
the  surface  of  a  small-scale  filament.  Numerical  simulations 
of  the  pulse  propagation  in  the  conditions  of  self-focusing 
and  normal  group- velocity  dispersion  [38]  also  propose  that 
the  short-pulse  supercontinuum  conical  emission  is  the  re¬ 
sult  of  four- wave  mixing  in  the  medium.  However,  according 
to  these  simulations,  conical  emission  should  also  exist  in 
the  Stokes  part  of  the  spectrum.  The  latter  fact  was  not  ob¬ 
served  in  the  experiments.  According  to  [2, 39]  the  conical 
emission  is  the  result  of  spatio-temporal  self-phase  modula¬ 
tion  of  a  femtosecond  pulse  propagating  in  the  conditions  of 
self-focusing  and  plasma  production.  The  temporal  phase  gra¬ 
dient  dcp/dz ,  which  defines  the  frequency  deviation  Sco  (r),  is, 
as  was  mentioned  earlier,  proportional  to  the  time  derivative 
of  the  electron  density:  Sco  (r)  =  d<p/ dr  ~  3 N&  (r)  / 3r.  At  the 
same  time,  the  spatial  phase  gradient,  dcp/dr ,  which  defines  the 
transverse  component  of  the  wave  vector  Skr ,  is  proportional 
to  the  radial  derivative  of  the  electron  density:  Skr  =  dcp/dr  ~ 

1 3 Ng  (r)  /dr\.  In  the  course  of  the  nonlinear  interaction  high 
spatio-temporal  gradients  of  free-electron  density  arise  sim¬ 
ultaneously  in  space  and  time.  As  a  result,  high-frequency 
spectral  components  propagate  at  a  larger  angle  to  the  propa¬ 
gation  direction. 

We  have  discussed  the  possible  origin  and  angular  dis¬ 
tribution  of  the  supercontinuum  blue  wing.  There  could  also 
exist  an  asymmetric  extension  to  the  Stokes  side  or  a  red 
shift  of  the  pulse  frequency  spectrum.  One  of  the  reasons  for 
the  Stokes  spectral  broadening  is  the  noninstantaneous  Kerr 
response  of  the  neutral  medium.  The  physical  mechanism 
of  it  is  a  nuclear  response  due  to  rotational  and  vibrational 
motions.  The  calculations  and  experiment  performed  in  [40] 
showed  that  the  response  time  of  the  rotational  Raman  mo¬ 
lecular  nonlinearity  in  molecular  gases  O2  and  N2  is  approxi¬ 
mately  70  fs.  The  combination  of  an  instantaneous  electronic 
Kerr  nonlinearity  and  a  delayed  rotational  Raman  molecu¬ 
lar  response  leads  to  a  steeper  leading  edge  in  the  positive 
time-dependent  change  of  the  refractive  index  as  compared 
with  the  trailing  edge.  As  a  result,  the  frequency  deviation 


Sco  (t)  =  —kzdAn/dz  is  larger  on  the  Stokes  side  and  the 
spectrum  exhibits  a  red-enhanced  part. 

In  the  conditions  of  a  strong  localization  of  energy  that 
enables  us  to  obtain  an  efficient  redistribution  of  input  laser 
pulse  energy  into  the  supercontinuum  energy,  there  exists  an 
additional  mechanism  that  enhances  the  blue  part  of  the  pulse 
spectrum  [41].  This  mechanism  is  temporal  self- steepening  of 
the  radiation  and  shock-wave  formation  at  the  trailing  edge 
of  the  pulse.  As  a  result,  the  positive  frequency  deviation 
Sco  (r)  =  —kzdAn/dz  at  the  back  of  the  pulse  is  larger  than 
the  absolute  value  of  the  negative  frequency  deviation  at  the 
front  of  the  pulse.  The  spectrum  acquires  a  large  blue  pedestal 
due  to  a  purely  positive  nonlinear  contribution  to  the  refrac¬ 
tive  index  from  the  Kerr  nonlinearity  [42-45].  The  effect 
is  the  stronger  the  shorter  is  the  input  pulse  duration.  Self- 
steepening  of  the  pulse  can  explain  the  formation  of  the  blue 
spectral  wing  in  the  propagation  of  ultra-short  pulses  through 
microstructured  fibers  [46].  In  these  optical  devices  a  steep 
trailing  edge  of  the  subpulse  formed  inside  the  fiber  may  origi¬ 
nate  due  to  the  joint  influence  of  higher-order  dispersion  terms 
(>  3)  and  the  Kerr  nonlinearity. 

In  summary,  the  effect  of  various  physical  mechanisms 
on  spectral  broadening  can  be  interpreted  in  terms  of  spatio- 
temporal  self-phase  modulation.  Indeed,  for  the  instantaneous 
electronic  Kerr  response,  the  ‘nonlinear’  phase  is  directly  pro¬ 
portional  to  the  intensity.  In  the  case  of  free-electron  gener¬ 
ation,  the  electron  density,  growing  as  a  function  of  the  in¬ 
tensity,  influences  the  ‘nonlinear’  phase  through  the  plasma 
contribution  to  the  refractive  index.  Stimulated  Raman  scat¬ 
tering  on  rotational  transitions  in  molecular  gases  results 
in  the  delayed  response  of  the  Kerr  nonlinearity  and  also 
leads  to  the  intensity-dependent  phase  of  the  pulse.  Self- 
steepening  increases  the  ‘nonlinear’  phase  gradients  through 
the  increase  of  intensity  gradients  on  the  pulse  fronts.  Such 
mechanisms  as  self-focusing  and  ionization  affect  not  only 
longitudinal  but  also  the  transverse  electric  field  distribution 
and,  therefore,  the  transverse  spatial  phase  of  the  pulse.  This 
gives  rise  to  such  phenomena  as  the  supercontinuum  conical 
emission. 

3  Experimental  setups  for  ultra-short  high-power 

pulse  generation 

In  the  experiments  we  used  a  Ti: sapphire  laser  sys¬ 
tem  (Fig.  la)  that  emits  three  simultaneous  beams.  The  system 
consists  of  an  oscillator  (Maitai,  28  fs  at  80  MHz,  Spectra 
Physics),  followed  by  a  stretcher  and  a  regenerative  amplifier 
(Spitfire,  Spectra  Physics),  which  works  at  1-kHz  repetition 
rate.  A  two-pass  amplifier  enhances  the  output  from  the  re¬ 
generative  amplifier.  The  amplified  beam  is  compressed  to 
40  fs  and  gives  a  maximum  energy  of  2  mJ/pulse  at  1  kHz.  On 
the  other  hand,  just  after  the  regenerative  amplifier,  a  pulse 
slicer  is  used  to  pick  up  one  pulse  every  100  ms.  This  10-Hz 
seed  beam  is  again  split  into  two  beams  by  a  beam  splitter. 
One  part  of  the  seed  beam  is  sent  to  a  two-pass  amplifier  and 
compressed  by  a  portable  compressor.  After  the  portable  com¬ 
pressor,  10-Hz  laser  pulses  with  an  energy  of  20mJ/pulse 
and  42-fs  duration  are  obtained.  The  second  10-Hz  seed  beam 
is  sent  to  a  four-pass  amplifier,  which  is  followed  by  a  vac¬ 
uum  compressor.  Through  a  10-m-long  pipe  with  a  1 -cm- thick 
CaF2  window  at  the  end,  our  vacuum  compressor  is  connected 
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4  Theoretical  model  for  the  nonlinear-optical  pulse 

transformation 

Transformation  of  the  shape  and  the  spectrum 
of  a  femtosecond  laser  pulse  results  from  the  nonlinear- 
optical  interaction  of  the  light  field  with  the  medium,  which 
is  defined  by  a  number  of  factors:  Kerr  self-focusing,  ion¬ 
ization  of  the  medium,  material  dispersion,  diffraction,  and 
self- steepening. 

Theoretical  models,  which  describe  the  effect  of  all  the 
factors  on  the  laser  pulse  propagating  in  transparent  di¬ 
electrics,  are  based  on  the  scalar  approximation  of  the  wave 
equation  for  the  light  field  q(x,y,  z,t)  in  the  dispersive 
medium  [47] : 


/  32  \  1  92  f  ,  ,  , 

W  +  a±)“-7^I 

—  00 

_  4jr  /  dj  32£>ni\ 

“  c2  V  3 1  3 12  )  ’ 


(1) 


where  pn\  is  the  nonlinear  polarization  of  the  neutral  atoms 
and  molecules  in  the  medium  and  j  is  the  current  of  the  free 
charges  produced  due  to  the  ionization  in  the  high-intensity 
laser  field. 

In  order  to  obtain  the  first-order  equation  relative  to  the 
propagation  coordinate  z,  let  us  represent  the  light  field 
q(x,y,  z,t)  in  the  form  of  the  wave  packet  with  the  carrier 
frequency  coq: 


a 


q(x,  y,  z,  t)  =  E(x,  y,  z,  t)  exp  i (coot  —  kz), 


(2) 


FIGURE  1  a  Ti: sapphire  laser  system  that  emits  three  beams  simultan¬ 
eously:  40  fs,  2mJ  per  pulse  at  1-kHz  repetition  rate,  42  fs,  20  mJ  per  pulse 
at  10-Hz  repetition  rate,  45  fs,  100  mJ  per  pulse  at  10-Hz  repetition  rate, 
b  Experimental  setup  used  for  supercontinuum  generation  in  water 

to  the  200-m-long  corridor.  At  the  exit  of  the  pipe,  we  can  have 
a  45 -fs  laser  pulse  and  up  to  100-mJ  pulse  energy.  The  typ¬ 
ical  spectrum  of  all  three  beams  is  centered  at  810  nm  with 
a  spectral  width  of  32  nm  (FWHM). 

For  the  experiments  on  the  supercontinuum  generation  in 
water  a  1-kHz  beam  with  a  maximum  energy  of  2mJ  per 
pulse  was  used  (Fig.  lb).  The  input  pulse  energy  was  con¬ 
trolled  by  a  half-wave  plate,  which  was  located  before  the 
compressor.  The  laser  pulse  was  focused  by  a  microscope 
objective  into  a  glass  water  cell.  The  transmitted  pulse  spec¬ 
trum  is  collimated  after  the  cell  by  an  achromatic  lens  and 
imaged  onto  the  slit  of  the  spectrometer.  Frequency  spec¬ 
tra  obtained  after  the  propagation  in  water  are  discussed 
in  Sect.  6. 


where  E(x,  y,  z,  t)  is  the  slowly  varying  complex  amplitude  of 
the  light  field.  Complex  amplitudes  of  the  nonlinear  polariza¬ 
tion  Pn i  and  the  field  current  J  at  the  frequency  qjq  are  given  by 
the  equations: 


Pnl  =  —Ank(I)E, 

IX 

dJ  e2  cno  d 

—  =  —JVc(I)E+-±-(a(I)E), 
at  rae  471  at 


(3) 

(4) 


where  no  is  the  unperturbed  refractive  index  of  the  medium, 
A /ik (7)  is  the  contribution  of  the  electronic  and  nuclear  re¬ 
sponse  of  the  Kerr  nonlinearity,  Ne(I)  is  the  free-electron 
density,  the  change  of  which  in  space  and  time  is  defined  by 
the  rate  equations  for  the  photo-ionization,  I  =  cno  \  E\ 2  /87T, 
and  mQ  and  e  are  the  electron  mass  and  charge,  respectively. 
The  absorption  coefficient  a (I)  that  describes  the  ionization 
energy  loss  is: 


a  —  l  ]mh(L>o 


(5) 


where  m  is  the  order  of  the  multiphoton  process. 

Following  [42],  we  can  obtain  from  (1)— (4)  the  nonlin¬ 
ear  equation  for  the  complex  amplitude  E(x,y,  z,t),  which 
in  the  retarded  coordinate  system  (r  =  t  —  z/ vg)  takes  the 
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form  [44, 45]: 
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The  first  term  on  the  right-hand  side  of  (6)  describes  the 
beam  diffraction;  the  second  and  the  third  terms  describe 
material  dispersion.  Note  that  in  the  conditions  of  strong 
self-focusing  and  self-steepening,  material  dispersion  should 
be  considered  up  to  the  third  order  or  higher,  because  the 
higher-order  dispersion  terms  influence  the  formation  of  sub¬ 
pulses  [45].  The  fourth  term  describes  the  contribution  of  the 
Kerr  nonlinearity  and  the  plasma  to  pulse  transformation.  The 
fifth  term  describes  the  pulse  energy  loss  due  to  the  ionization. 

Equation  (6)  for  the  complex  amplitude  of  the  electric 
field  E ,  unlike  the  wave  equation  in  the  slowly  varying 
envelope  approximation  (SVEA)  [41],  contains  the  opera¬ 
tor  [1  zb  (i/coo) (d/ dr)],  which  appears  after  the  substitution 
of  (2)-(4)  into  (1)  and  retaining  the  higher-order  terms  in 
comparison  with  the  SVEA.  This  approximation  with  the 
higher-order  correction  terms  is  called  the  slowly  evolving 
wave  approximation  (SEWA)  [42].  It  was  shown  in  [42] 
that  in  the  framework  of  the  SEWA  the  equation  for  the 
complex  amplitude  of  the  electric  field  E  accurately  de¬ 
scribes  light-pulse  propagation  down  to  the  single-cycle 
regime.  For  this  pulse  the  spectral  width  is  comparable  to 
the  carrier  frequency  coo.  When  applied  to  the  problem  of 
nonlinear-optical  transformation  of  an  ultra-short  high-power 
pulse,  the  SEWA  adequately  reproduces  the  formation  of 
a  wide  broadened  spectrum  of  the  radiation.  In  terms  of 
mechanisms  that  contribute  to  the  modification  of  ultra- 
short  radiation,  the  SEWA  describes  self-steepening  of  the 
pulses  and  subpulses  arising  in  the  course  of  propagation 
as  well  as  the  shock-wave  formation  at  the  back  of  the 
pulse. 

The  Kerr  nonlinearity  of  neutrals  far  from  the  atomic  and 
molecular  resonances  is  defined  by  an  anharmonic  response 
of  the  bound  electrons  and  the  stimulated  Raman  scattering 
on  rotational  transitions  of  molecules.  In  the  case  of  fem¬ 
tosecond  pulses,  the  electronic  response  may  be  considered 
as  instantaneous.  The  contribution  from  the  stimulated  Raman 
scattering  is  nonstationary,  as  shown  theoretically  in  [48]  and 
experimentally  and  theoretically  in  [40].  In  [16],  the  response 
function  H{t)  measured  in  [40]  was  approximated  based  on 
the  damped  oscillator  model  by  the  following  equation: 

Hit)  =  0{t)£22  exp(— 71/2)  sin  (/If)  /A  ,  (7) 


where  0{t)  is  the  Heaviside  function  and  A  —  v/ Q-  —  A1  /4. 
The  fitting  parameters  have  the  following  values:  Q  — 
20.6  THz  and  T  =  26THz.  The  corresponding  characteris¬ 
tic  times  Ti  =  2 /T  and  r2  =  1  /A  are  X\  —  11  fs  and  t2  = 
62.5  fs.  These  values  are  comparable  with  the  laser-pulse 
duration  used  in  the  recent  laboratory  experiments  on  the 
observation  of  filamentation  and  supercontinuum  genera¬ 
tion  [9-14, 18-20]. 


Finally,  the  nonlinear  contribution  to  the  refractive  index 
from  the  neutral  molecules  is  considered  in  the  form: 

t 

Ank(t)  =  (l-g)^n2\E\2  +  g  J  l-n2\E{r,t')\2H{t-t')At'. 

—  00 

(8) 


The  nonlinear  contribution  to  the  refractive  index  A nv 
(a,  y,  z,  t)  from  the  laser-produced  plasma  is  given  by: 


A nv(x,  y,  z,  t) 


a)gp,  [ 0  /  vA 

2no(a>l  +  v2c)  \  coqJ 


(9) 


where  cov(x,  y,  z,  t)  =  e2N&(x,  y,  z,  t)/mG  is  the 

plasma  frequency,  vc  =  Nq  veac  is  the  effective  electron  colli¬ 
sion  frequency  with  the  root-mean-square  electron  velocity  ve 
and  the  electron  collision  cross  section  crc,  and  No  is  the  dens¬ 
ity  of  neutrals. 

The  free-electron  density  Ne(x,  y,  z,  t)  depends  on  the 
spatial  coordinates  and  time  according  to  the  kinetic  equation: 

=  fl(|£|2)(Af0  -  Ne)  +  V,- A/e  -  /we2,  (10) 

ot 

where  No  and  R(|£|2)  are  the  density  and  the  ionization 
rate  of  neutrals,  respectively  and  is  the  radiative  electron- 
recombination  coefficient.  The  avalanche  ionization  fre¬ 
quency  is: 


1  e2£2 

W  2me(col  +  v2) 


(11) 


where  W  is  the  ionization  potential  in  the  case  of  a  gaseous 
medium  or  the  band  gap  in  the  case  of  condensed  media. 

Equation  (10)  is  written  in  the  single-ionization  approx¬ 
imation.  This  approximation  can  be  applied  to  the  study  of 
ultra-short  pulse  propagation  in  air  since  the  ionization  de¬ 
gree  Nq / No  is  less  than  0.001  [49, 50].  In  water  the  ionization 
degree  reaches  0.02  due  to  the  essential  contribution  of  the 
avalanche  ionization  even  in  the  case  of  an  ultra-short  45-fs 
pulse  [51]. 

In  the  case  of  a  femtosecond  pulse  propagation  in  atmo¬ 
spheric  density  gases,  the  effective  electron-neutral  collision 
frequency  vc  is  much  smaller  than  the  laser  frequency  coo  [52], 
and  (9)  for  the  plasma  contribution  to  the  refractive  index 
takes  the  form: 


A np(x,  y,  z,  t) 


ft)2(x,  y,  z,  t) 

2n0a> l 


(12) 


In  the  atmospheric  air  the  large  amount  of  free  electrons 
comes  from  the  oxygen  molecules,  the  ionization  potential  of 
which  Wo2  =  12.1  eV  is  smaller  than  the  ionization  potential 
of  the  nitrogen  molecules  W^2  =  15.6  eV  [53].  To  calculate 
the  ionization  rate  R  (|  E\2)  in  (10)  we  used  the  model  [54]  for 
the  ionization  of  a  hydrogen-like  atom  in  a  linearly  polarized 
electric  field  E: 


R  =  I  Crc*/*  I2  flmEjJ  ~ 
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Here  y  =  cooV2W  /E  is  the  Keldysh  parameter,  n*  =  Z/ *J2Ei 
is  the  effective  main  quantum  number,  ~  n*  —  1  is  the  ef¬ 
fective  orbital  quantum  number,  Z  is  the  effective  ion  charge, 
/  is  the  orbital  quantum  number,  and  m  is  the  projection 
of  the  orbital  quantum  number  on  the  direction  of  the  elec¬ 
tric  field  vector.  The  coefficients  Cn*i*  and  //m  as  well  as 
the  functions  Am  (coo,  y)  and  g(y)  can  be  found  in  [55]. 
The  values  of  effective  charges  for  the  ionization  of  oxy¬ 
gen  and  nitrogen  molecules  are  calculated  in  [55]  by  fitting 
(13)  to  the  experimental  data  on  molecular  ion  yields.  They 
are  equal  0.53  and  0.9  for  oxygen  and  nitrogen  molecules, 
respectively. 

The  pulse  shape  and  beam  profile  at  the  femtosecond  laser 
system  output  are  close  to  Gaussian;  therefore,  the  initial  dis¬ 
tribution  of  the  electric  field  complex  amplitude  in  the  re¬ 
tarded  coordinate  system  is  described  by: 


E(x,  y,  z  =  0,  t)  =  E0 


_T0_ 

Tp(5) 


xexpl 


x2  +  y2 
2a20 


r2  .  kr2  ,5t2N 

— ~ - b  i - b  i - 

2r2(<5)  2  Rf  2 


(14) 


where  ao  is  the  input  beam  pulse  radius  at  1/e  intensity  level, 
Rf  is  the  geometrical  focusing  distance,  and  To  is  a  half-pulse 
duration  of  a  transform-limited  pulse  at  1  /  e  intensity  level. 
The  parameter  To  is  related  to  the  full  width  at  half  maximum 
pulse  duration  tfwhm,  which  is  usually  measured  in  the  ex¬ 
periments,  as  To  =  TpwHM/(2\/ln  2).  If  at  the  output  of  the 
compressor  the  pulse  is  chirped,  then  the  initial  pulse  dura¬ 
tion  Tp(<$)  depends  on  the  parameter  8(zp(8  =  0)  =  To),  which 
characterizes  the  initial  phase  modulation  of  the  pulse  with 
a  constant  spectral  width  [41, 56]. 

The  solution  to  the  self-consistent  problem  (7)-(10),  (13) 
with  the  initial  conditions  (14)  defines  the  complex  amplitude 
of  the  electric  field  E(x,  y,  z,  t)  after  the  self-transformation 
of  the  pulse  in  the  nonlinear  dispersive  medium.  The  fre¬ 
quency  angular  spectrum  S(0X,  0y,  A co,  z)  at  a  distance  z  is 
calculated  through  the  spatio-temporal  Fourier  transform  as: 

5  (0X,  9y,  Aw,  z )  =  I  U(0X,  9y,  Aw,  z) | 2  , 

U(dx,0y,  Aw,  z)=J  dxdye-Wxkxc~i9yky 

x  J  dre~'Aa>rE(x,y,z,r),  (15) 

where  0X  =  kx/k,  0y  =  ky/k  are  the  angles  at  which  differ¬ 
ent  frequency  components  A  co  propagate  in  the  medium.  In 
the  visible  part  of  the  frequency  spectrum  and  for  the  radially 
symmetric  pulse  distribution,  the  dependence  of  the  angle  0  = 
0X  =  0y  on  wavelength  defines  the  supercontinuum  cone  emis¬ 
sion.  In  order  to  represent  the  conical  emission  in  the  coor¬ 
dinates  6(k),  the  following  transition  from  spectral  frequency 
to  spectral  wavelength  is  made:  k  =  2nc/  (coo  +  A  co),  where 
coo  is  the  laser  central  frequency.  For  the  spatio-temporal  fre¬ 
quency  deviation  as  the  function  of  time  and  radial  coordinate, 
the  corresponding  spectral  wavelengths  are  calculated  accord¬ 
ing  to  the  expression: 


k  =  2jxc/  (coq  +  3 <p/ 3t)  , 


(16) 


where  cp  is  the  phase  of  the  complex  amplitude  of  the  electric 
field  E(r,r). 

5  Supercontinuum  generation  in  air 

Supercontinuum  generation  in  air  is  closely  con¬ 
nected  to  the  phenomenon  of  high-power  femtosecond  laser 
pulse  filamentation.  This  phenomenon  is  characterized  by  the 
stabilization  of  the  radiation  parameters,  which  are  available 
for  measurements.  Among  such  parameters  are  the  filament 
diameter  and  energy.  Both  values  are  defined  from  the  nar¬ 
row  near-axis  part  of  the  fluence  distribution.  In  addition,  the 
peak  intensity  and  the  maximum  plasma  density  in  the  fila¬ 
ment  preserve  the  same  values  of  ~  4  x  1013  W/cm2  [53,57] 
and  1016  cm-3  [50],  respectively.  The  clamping  of  the  peak  in- 
tensity  inside  the  filament  to  a  maximum  value  can  be  readily 
understood  [53]  to  be  due  to  a  balance  between  the  Kerr  self- 
focusing  and  the  plasma  defocusing  effects.  It  has  been  shown 
in  gases  by  observing  the  plateau  in  high-harmonic  generation 
for  Xe  [58]  and  by  observing  fluorescence  spectra  emitted 
from  excited  nitrogen  molecules  and  molecular  ions  [57]. 

Fluence  is  an  integrated  characteristic  of  the  propagation 
and  its  distribution  does  not  reveal  the  mechanism  of  super¬ 
continuum  generation.  In  spite  of  the  fact  that  the  fluence  is 
stabilized,  a  spatio-temporal  intensity  distribution  in  the  fil¬ 
ament  is  strongly  nonstationary  and  experiences  significant 
changes  along  the  propagation.  The  narrow  near-axis  part 
(100  |xm)  of  the  pulse  is  fed  from  the  wide  (of  the  order  of 
or  larger  than  the  input  beam  diameter)  background  of  the 
radiation,  which  contains  around  90%  of  the  total  pulse  en¬ 
ergy.  The  dynamic  process  of  the  energy  replenishment  in  the 
filament  in  the  conditions  of  the  Kerr  nonlinearity,  the  laser- 
produced  plasma,  natural  beam  divergence,  material  disper¬ 
sion,  and  self- steepening  leads  to  dramatic  transformation  of 
the  intensity  and  phase  of  the  light  field  on  the  whole  trans¬ 
verse  beam  section,  which  exceeds  the  input  beam  diameter. 
This  results  in  the  superbroadening  of  the  frequency  spectrum 
and  conical  emission  generation.  In  order  to  adequately  model 
the  self-transformation  of  the  pulse  shape  and  spectrum  as 
well  as  to  obtain  quantitative  agreement  with  the  experimental 
data  on  the  amount  of  energy  in  the  filament,  the  peak  electron 
density,  the  width  of  the  frequency  spectrum,  and  the  conical 
emission  angles,  both  the  small-scale  changes  in  the  filament 
center  and  the  large-scale  changes  in  the  background  should 
be  taken  into  account  in  the  simulations. 

The  experiment  on  the  study  of  supercontinuum  cone 
emission  was  performed  using  a  Ti: sapphire  chirped  pulse 
amplification  laser  system  generating  250-fs  (FWHM)  laser 
pulses  with  the  central  wavelength  of  800  nm  (the  geometry 
of  the  setup  was  similar  to  that  shown  in  Fig.  la  except  that 
the  smaller  bandwidth  from  the  oscillator  and  the  stretcher  as 
well  as  gain  narrowing  allowed  us  to  generate  only  the  pulses 
longer  than  lOOfs  FWHM).  The  beam  radius  at  the  output 
of  the  grating  compressors  was  ao  =  3.5  mm  (1/e  intensity 
level).  The  collimated  beam  was  sent  into  the  long  hallway  by 
means  of  several  transport  mirrors  and  propagated  freely  to 
the  end  of  the  hallway  for  a  total  distance  of  1 1 1  m.  In  order 
to  study  the  spatial  structure  of  the  white  light  accompanying 
filamentation,  the  cone  angle  of  the  various  colors  in  the  coni¬ 
cal  emission  was  measured.  In  the  experimental  setup  we  used 
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a  4-mm-diameter  aperture  that  was  centered  on  the  beam  axis 
and  a  razor  blade  that  blocked  off  the  filament.  The  conical 
emission  was  measured  on  a  white  screen  placed  at  20  m  after 
the  aperture.  An  interference  filter  with  10-nm  bandwidth  was 
placed  after  the  razor  blade.  The  radii  of  the  colored  half-rings 
were  measured  on  the  screen  at  the  following  wavelengths: 
500,  550,  600,  650,  700,  and  750  nm.  No  ring  was  observed 
at  X  >  800  nm.  The  half-angle  Ox  of  the  conical  emission  was 
defined  as  the  ratio  of  the  radius  r  of  the  ring  to  the  distance 
D  from  the  aperture.  Experimental  results  of  the  conical  emis¬ 
sion  angles  as  a  function  of  wavelength  are  shown  in  Fig.  2. 
The  pulse  energy  for  this  measurement  was  10  mJ. 

To  simulate  the  propagation  we  used  the  system  of  equa¬ 
tions  (7)-(8),  (10),  and  (12)  with  an  input  spatio-temporal 
distribution  of  the  electric  field  in  the  form  of  (14),  where 
r0  =  150fs  corresponds  to  250fsFWHM,  intensity  7o  = 
10nW/cm2,  beam  radius  <20  =  3.5  mm,  diffraction  length 
Zd  =  ka^  =  96  m,  and  <5  =  0.  The  pulse  energy  was  10  mJ  and 
the  peak  power  Fpeak  =  38  GW  =  6.3 Pcr,  where  Pcr  is  the 
critical  power  for  self-focusing  in  air.  The  calculations  were 
performed  on  a  grid  with  nonequidistant  steps  in  the  radial  di¬ 
rection  [59].  The  number  of  grid  steps  was  2048  in  the  time 
domain,  7000  along  the  propagation  coordinate  z,  and  512  in 
the  radial  direction. 

Redistribution  of  the  intensity  in  the  transverse  spatial  di¬ 
rection  and  in  time  is  shown  in  Fig.  3,  where  the  scale  in 
the  plane  (r,  t)  is  chosen  so  that  the  equal-intensity  contours 
at  the  beginning  of  propagation  (z  =  0)  are  concentric  rings 
(Fig.  3a).  At  the  start  of  the  filament  (z  =  0.29zd  ~  28  m)  the 
peak  intensity  of  the  pulse  is  attained  in  the  temporal  slice 
with  t  >  0  due  to  the  delayed  response  of  the  Kerr  nonlin¬ 
earity  (8)  (Fig.  3b).  The  value  of  the  peak  intensity  reaches 
~4x  1013  W/cm2  and  the  ionization  starts  to  defocus  the 
trailing  part  of  the  pulse.  At  z  —  0.32zd  ~  31  m  (Fig.  3c)  the 
intensity  maximum  shifts  towards  the  leading  part  of  the 
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FIGURE  2  Dependence  of  the  conical  emission  angle  on  wavelength  in 
the  blue  wing  of  the  pulse  frequency  spectrum.  Experimental  data  [38]  are 
shown  by  the  black  dots  and  the  simulation  results  by  the  solid  curve.  Both  in 
the  experiment  and  in  the  simulations  the  pulse  duration  to  =  150  fs  (250  fs 
FWHM),  the  pulse  energy  is  lOmJ  (the  ratio  of  the  pulse  peak  power  to 
the  critical  power  for  self-focusing  Ppeak /  Pa  =  6.3),  the  input  beam  radius 
ao  =  3.5  mm,  and  the  propagation  distance  z  =  40  m 


FIGURE  3  Spatio-temporal  distribution  of  the  laser  pulse  intensity  7(r,  r) 
at  different  distances  z  from  the  laser  system  output.  The  equal-intensity 
contours  are  plotted  for  the  following  intensity  values:  In  =  Iq2'\  where 
n  =  - 3,  -2, ...  ,  +9,  7o  «  1011  W/cm2.  a  z  =  0;  b  z  =  0.29 zd  «  28  m; 
c  z  =  0.32zd  ~  31  m;dz  =  0.4zd  «  39  m 


pulse;  its  value  reaches  6  x  1013  W/cm2.  At  the  trailing  part 
of  the  pulse  rings  are  formed.  By  z  =  0.4zd  ~  39  m  (Fig.  3d) 
the  diameter  of  the  outer  ring  exceeds  the  input  beam  diam¬ 
eter  by  a  factor  of  two.  At  the  leading  part  we  can  see  the  very 
intense  slices  between  z  ~  —  80  fs  and  r  ~  —  50  fs  located 
in  the  narrow  near-axis  part  within  the  transverse  diameter 
of  less  than  300  pun.  These  slices  form  the  ‘filament’  itself, 
i.e.  the  structure  that  appears  to  a  human  eye  as  a  string  of 
light  foci.  Formation  of  a  dynamic  multipeak  structure  in  the 
spatio-temporal  domain  of  the  pulse  is  revealed  as  the  spectral 
broadening  in  the  frequency  angular  domain. 

Because  of  the  correct  consideration  of  the  low-intensity 
background  of  the  radiation,  we  were  able  to  perform  a  quan¬ 
titative  comparison  between  the  conical  emission  angles  ob¬ 
tained  in  the  experiment  (Fig.  2,  symbols)  and  in  the  simula¬ 
tions  (Fig.  2,  solid  curve).  Note  that  in  the  previous  theoretical 
studies  [2, 39]  only  the  relative  angles  Ox/Oo ,  where  0q  is  the 
input  beam  divergence,  could  be  related  to  the  experimen¬ 
tal  data.  The  quantitative  discrepancy  between  the  simulated 
and  experimentally  observed  conical  emission  angles,  which 
reaches  ~  20%  for  500-nm  wavelength,  is  due  to  the  fact  that 
for  the  numerical  study  in  these  particular  experimental  con¬ 
ditions  we  were  not  able  to  include  material  dispersion  in  the 
simulations.  At  the  same  time  it  has  been  shown  by  Golubtsov 
et  al.  [39]  that  the  inclusion  of  the  material  dispersion  leads  to 
the  increase  of  the  conical  emission  angles. 

For  the  understanding  of  the  spectral  superbroadening 
mechanism  we  study  the  sources  of  the  supercontinuum  in 
the  pulse.  In  Fig.  4  we  present  the  map  of  spectral  wave¬ 
lengths  corresponding  to  the  frequency  deviation  800  (r,  r)  = 
dcp(r,  z)  /dr  (see  (16)).  The  distribution  of  wavelengths  is 
presented  together  with  the  spatio-temporal  intensity  distribu¬ 
tion  7(r,  r)  (equal-intensity  contours  are  shown  by  solid  lines). 
The  upper  inset  in  Fig.  4  shows  the  distribution  of  the  elec- 
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FIGURE  4  Spatio-temporal  distribution 
of  the  supercontinuum  sources.  The  solid 
lines  correspond  to  the  equal-intensity 
contours  plotted  for  the  following  in¬ 
tensity  values:  In  =  Io2n,  where  n  = 
-3,  —2, . . .  ,  +9,  70  «  1013  W/cm2.  The 
upper  inset  corresponds  to  the  tempo¬ 
ral  distribution  of  the  pulse  intensity 
{solid  line )  and  phase  (< dashed  line )  at 
r  =  13  |xm;  the  left  and  the  right  in¬ 
sets  correspond  to  the  radial  distribu¬ 
tion  of  the  intensity  {solid  line )  and  the 
phase  {dashed  line)  at  r  =  —1 40  fs  and 
r  =  217  fs,  respectively.  Propagation  dis¬ 
tance  z  —  0.76zd  =  73  m.  Pulse  param¬ 
eters  are  the  same  as  in  Fig.  3.  Radial 
position  r  =  13  pan  and  temporal  pos¬ 
itions  r  =  — 140  fs  and  r  =  217fs  are 
indicated  by  the  white  dashed  lines  on 
the  map 
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trie  field  intensity  (solid  line)  and  the  phase  (dashed  line)  as 
a  function  of  time  for  the  radial  position  r  =  73  |xm  indicated 
by  the  white  dashed  horizontal  line  on  the  map.  The  left  and 
the  right  insets  in  Fig.  4  show  the  distribution  of  the  electric 
field  intensity  (solid  line)  and  the  phase  (dashed  line)  as  the 
function  of  the  radial  position  at  the  leading  (t  =  — 140  fs)  and 
the  trailing  (r  =  217  fs)  edges  of  the  pulse.  The  temporal  pos¬ 
itions  are  indicated  by  the  white  dashed  vertical  lines  on  the 
map. 

In  Fig.  4  only  the  narrow  near-axis  part  of  the  pulse  distri¬ 
bution  with  the  diameter  of  the  order  of  300  |xm  is  presented. 
At  this  propagation  distance  z  =  0.76zd  ~  73  m  the  intensity 
distribution  reveals  two  peaks.  The  peak  that  is  formed  in 
the  leading  edge  of  the  pulse  r  =  —  140fs  is  due  to  the  self- 
focusing  in  air  and  the  peak  at  the  trailing  part  r  =  180  fs 
is  formed  due  to  the  joint  contribution  of  the  refocusing  and 
the  delayed  Kerr  response.  Two  distinct  ring  sources  of  the 
short- wavelength  components  in  the  frequency  spectrum  can 
be  seen  on  the  map.  Each  ring  source  is  represented  by  two 
white  points  located  symmetrically  relatively  to  the  beam  axis 
r  =  0.  The  first  is  at  the  leading  edge  of  the  pulse,  r  =  — 140  fs, 
and  its  radius  is  73  |xm.  This  is  the  region  where  the  Kerr- 
induced  convergence  is  replaced  by  the  plasma-induced  diver¬ 
gence  of  the  high-intensity  slices.  The  time  dependence  of  the 
phase  along  the  white  dashed  line  r  =  73  |xm  shows  the  strong 
jump  slightly  before  r  =  — 150  fs  and,  hence,  the  positive  fre¬ 
quency  deviation  8co  >  0  (upper  inset  in  Fig.  4).  Note  that,  at 
the  same  moment  r  =  —  150fs,  the  radial  phase  distribution 
reveals  a  large  gradient  dip /dr  at  the  ring  r  =  73  pun  (left  in¬ 
set  in  Fig.  4),  which  means  that  the  new  high  frequencies  born 
on  this  ring  diverge  at  large  angles  to  the  propagation  axis. 
Actually,  this  is  how  the  conical  emission  originates  from 
the  nonlinear  transformation  of  the  pulse  in  the  conditions  of 
self-focusing  and  the  plasma  defocusing.  The  second  super¬ 
continuum  source  starts  from  the  ring  at  the  trailing  edge  of 
the  pulse  r  =  180  fs,  r  =  68  pm  and  persists  over  the  whole 
back  front  r  >  180  fs.  Here  the  positive  frequency  deviation 
8co  >  0  is  due  to  the  shock  formation  at  the  back  of  the  pulse, 
which  can  also  be  seen  in  the  upper  inset  of  Fig.  4,  where  the 


phase  derivative  dip /dr  is  positive  at  r  >  187  fs.  In  the  radial 
phase  distribution  at  the  back  of  the  pulse  (r  =  217  fs,  see  the 
right  inset  in  Fig.  4)  the  absolute  value  of  the  phase  gradient 
dip /dr  is  much  smaller  than  in  the  front  of  the  pulse.  In  add¬ 
ition,  there  are  no  intensity  rings  corresponding  to  the  radial 
phase  jumps  (compare  the  radial  intensity  distributions  shown 
by  the  dashed  lines  in  the  left  and  the  right  insets  in  Fig.  4). 
Hence,  the  spatial  distribution  of  the  high-frequency  compo¬ 
nents  born  due  to  self- steepening  at  the  trailing  edge  of  the 
pulse  is  not  in  the  form  of  conical  rings. 

Figure  5  clearly  demonstrates  the  role  of  self-steepening 
in  the  supercontinuum  generation.  In  Fig.  5  the  on-axis  fre¬ 
quency  spectrum  for  the  case  without  self- steepening,  i.e. 
without  the  operator  [1  =b  (i /coo)(d/ dr)]  in  (7),  is  plotted  by  the 
solid  line  and  the  spectrum  obtained  with  the  inclusion  of  this 
operator  is  plotted  by  the  dotted  line.  For  the  on-axis  spectrum 


FIGURE  5  Frequency  spectrum  calculated  on  the  beam  axis  without  con¬ 
sideration  of  the  self-steepening  term  ±(//<u o)(3/3r)  in  (6)  {solid  line)  and 
with  the  inclusion  of  this  term  so  that  (u>oto)-1  =  0.007  {dotted  line).  In  the 
simulations  the  pulse  duration  was  tq  =  90  fs  (150fs  FWHM),  the  ratio  of 
the  peak  power  to  the  critical  power  for  self-focusing  Tpeak /  Ar  =  6,  the  in¬ 
put  beam  radius  ao  =  0. 18  mm,  and  z  =  Zd-  The  intensity  is  normalized  to  the 
peak  input  spectral  intensity 
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FIGURE  6  Frequency  angular  spec¬ 
trum  (a,  b)  and  spatio-temporal  distribu¬ 
tion  (c,  d)  of  the  pulse  at  two  propagation 
distances:  a,  c  z  =  0.43zd;  b,  d  z  =  0.6zd- 
In  the  simulations  the  pulse  duration  was 
r0  =  150  fs  (250  fs  FWHM),  the  ratio  of 
the  peak  power  to  the  critical  power  for 
self-focusing  /peak/ Pa  =  6.3,  and  the  in¬ 
put  beam  radius  ao  =  0. 18  mm 


at  this  propagation  distance  the  effect  of  self- steepening  leads 
to  the  increase  of  the  blue  wing  by  250  nm  at  10-4  level  of  the 
maximum  intensity. 

If  the  frequency  spectrum  is  integrated  over  the  whole 
transverse  aperture,  the  relative  contribution  of  self- steepening 
to  the  spectral  broadening  decreases  [45].  The  less  pro¬ 
nounced  effect  of  self- steepening  in  the  case  of  the  integrated 
spectrum  is  associated  with  the  contribution  of  the  high- 
frequency  spectral  components,  which  are  generated  at  the 
off-axis  positions  mainly  due  to  the  laser-induced  plasma 
production.  This  is  confirmed  by  the  frequency  angular  spec¬ 
trum  shown  in  Fig.  6a  and  b  for  the  two  distances  z  =  0.43zd 
and  z  =  0.6zd  along  the  filament.  The  corresponding  spatio- 
temporal  intensity  distribution  are  shown  in  Fig.  6  (c,d),  re¬ 
spectively.  The  same  figure  demonstrates  spatio-temporal 
intensity  distributions  at  the  corresponding  propagation  dis¬ 
tances.  At  the  beginning  of  the  filamentation  the  blue  wing 
is  mainly  formed  in  the  ring  at  the  leading  edge  of  the  pulse 
(compare  with  Fig.  4).  Therefore,  the  conical  divergence  of 
high-frequency  components  dominates  in  the  pulse  spectrum 
(Fig.  6a).  Later  on,  the  pulse  breaks  up  into  the  two  subpulses. 
At  the  trailing  edge  of  the  second  subpulse  (t  =  140  fs)  the 
self- steepening  is  pronounced  on  the  beam  axis  (Fig.  6d).  As 
a  result,  the  high-frequency  components  that  appear  in  the 
pulse  spectrum  do  not  have  the  angular  divergence  and  propa¬ 
gate  mainly  along  the  axis  (Fig.  6b). 

Thus,  the  supercontinuum  generation  in  air  strongly  de¬ 
pends  on  the  spatio-temporal  phase  gradients  of  the  complex 
amplitude  of  the  electric  field.  In  Sect.  8  we  will  show  how  the 
phenomenon  of  long-range  filamentation  in  air  can  allow  us  to 
control  this  spatio-temporal  gradient  and  provide  the  most  ef¬ 
fective  redistribution  of  energy  to  the  supercontinuum  energy 
in  the  prescribed  position  in  the  propagation. 

6  Supercontinuum  generation  in  water 

Supercontinuum  generation  in  water  occurs  al¬ 
ready  at  low  (megawatt)  peak  input  powers,  since  the  nonlin¬ 
ear  refractive  index  is  larger  than  that  of  air.  The  first  set  of 
experiments  in  water  was  performed  with  the  pulses  generated 


by  the  laser  system  similar  to  the  one  shown  in  Fig.  1  a  but  with 
narrower  bandwidth  components  that  allowed  us  to  generate 
only  the  pulses  longer  than  100  fs  FWHM.  The  pulse  duration 
was  170  fs  FWHM  and  the  energy  varied  from  0.2  to  200  |xJ. 
The  beam  diameter  at  1  /  e2  level  was  4  mm.  According  to  the 
scheme  in  Fig.  lb,  the  beam  was  focused  with  a  119-cm  lens 
into  a  cell  filled  with  water.  Further  details  of  the  experimental 
setup  can  be  found  in  [60].  Spectra  obtained  at  low  and  high 
input  laser  energies  are  shown  in  panels  (a)  and  (b)  of  Fig.  7 
respectively.  Each  spectrum  is  normalized  to  its  maximum  at 
the  central  wavelength  of  800  nm.  While  the  spectra  at  low 
energies  show  a  symmetrical  broadening  towards  lower  and 
higher  frequencies,  we  observe  the  typical  long  pedestal  on 
the  blue  side  at  higher  pulse  energies  (note  that  the  bandwidth 
of  the  spectrometer  limited  the  observations  on  the  red  side 
to  wavelengths  below  1015  nm).  The  minimum  wavelength  of 
the  blue  spectral  wing  remains  almost  constant  for  an  increase 
in  the  pulse  energy  by  two  orders  of  magnitude  from  0.4  to 
200  |xJ  (this  wavelength  is  well  within  the  bandwidth  of  the 
blue  filter  used  in  the  experimental  setup).  In  order  to  illustrate 
this  more  clearly,  we  present  in  Fig.  8  the  maximum  positive 
frequency  shift  as  a  function  of  the  average  input  power.  Simi¬ 
lar  observations  of  a  cut-off  wavelength  on  the  blue  side  have 
been  made  in  the  supercontinuum  spectra  of  chloroform  and 
glass  [60]. 

For  a  qualitative  understanding  of  these  observations  it  is 
sufficient  to  analyze  which  parameters  of  the  laser  pulse  are 
mainly  responsible  for  the  frequency  shifts  in  the  supercontin¬ 
uum  generation.  Considering  plane  waves  propagating  in  the 
nonlinear  medium,  the  frequency  deviation  at  a  certain  propa¬ 
gation  distance  z  can  be  written  as  [60] : 


8(0  (t)  =  —kz 


3  (Ant  +  Aflp) 
3r 


(17) 


where  An^  and  A nv  are  the  contributions  to  the  refractive  in¬ 
dex  from  the  Kerr  nonlinearity  and  the  plasma,  given  by  (8) 
and  (9),  respectively.  For  a  simple  estimate  we  can  consider 
the  Kerr  nonlinear  response  of  water  as  instantaneous  (g  =  0 
in  (8))  and  take  into  account  that  the  ionization  of  water  starts 
from  the  multiphoton  process  of  the  order  of  m  =  5  for  the 
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FIGURE  7  Spectra  obtained  for  the  propagation  of  a  Ti: sapphire  laser  pulse 
with  the  duration  170fs  FWHM  in  water  at  various  input  pulse  energies 
(see  legend):  a  at  low  laser  energies  (0.4- 1.8  pi);  b  at  higher  laser  energies 
(2-200  |xJ) 


pulse  with  the  central  wavelength  of  800  nm.  Then  (17)  can  be 
rewritten  in  the  form  (e.g.  [62]): 

8co{x)  =  -a^-^-+bIm  (x)  ,  (18) 

dr 

where  a  =  Azzn^kzlc,  b  =  2n  e2 /3(m) Nokz/ (rae&>o)  are  ^de¬ 
pendent  of  time  and  /?(m)  is  the  ra-photon  absorption  coeffi¬ 
cient  of  water. 

At  low  laser  energies  (0.4- 1.8  pJ)  the  symmetrical  spec¬ 
tral  broadening  is  due  to  the  first  term  on  the  right-hand  side 
of  (18)  (Fig.  7a).  At  higher  laser  energies  the  second  term  on 
the  right-hand  side  of  (18)  starts  to  contribute  to  the  spectral 
broadening  and  creates  the  asymmetric  blue  wing.  From  2  to 
200  pJ  the  ratio  of  the  pulse  peak  power  Fpea k  to  the  critical 
power  for  self-focusing  in  water  PCT  =  4.4  MW  [61]  ranges 
from  Ppeak/ Per  =  2.6  to  Fpeak/Fcr  =  260.  The  intensity  growth 
due  to  self-focusing  is  stopped  by  the  ionization.  The  thresh¬ 
old  intensity,  at  which  the  ionization  stops  self-focusing,  is 
a  function  of  the  medium  and  above  a  certain  pulse  energy 
nearly  does  not  depend  on  the  specifics  of  the  pulse  propa¬ 
gation  [53].  Thus,  the  constant  maximum  positive  shift,  ob¬ 
served  in  Figs.  7b  and  8,  is  the  experimental  evidence  of  the 
existence  of  such  a  threshold  intensity  or,  in  other  words,  of 
the  peak-intensity  clamping. 

In  (17)  and  (18)  we  have  specified  the  main  contributions 
to  the  frequency  deviation  8co  in  the  bulk  medium.  How¬ 
ever,  since  this  deviation  is  a  function  of  the  intensity  gra- 
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FIGURE  8  The  maximum  positive  frequency  shift  A<umax  as  a  function  of 
the  average  input  power,  obtained  for  the  propagation  of  a  170-fs  FWHM 
pulse  in  water.  The  average  power  of  1  mW  corresponds  to  the  energy  of  1  |iJ 
per  pulse 


dient,  all  processes  in  the  medium  that  affect  the  intensity 
gradients  simultaneously  affect  the  distribution  of  frequen¬ 
cies  in  the  supercontinuum.  One  of  such  factors  is  material 
dispersion.  Kolesik  et  al.  [63]  compared  the  extension  of  the 
supercontinuum  wings  generated  in  air  in  the  conditions  of 
material  dispersion  approximated  up  to  the  second  order  and 
the  third  order  with  the  extension  generated  when  full  chro¬ 
matic  dispersion  of  air  is  taken  into  account.  The  error  in¬ 
troduced  by  the  incomplete  description  of  the  dispersion  is 
revealed  at  the  level  of  10-4  of  the  power  spectrum  amplitude 
at  the  blue  side.  Thus,  material  dispersion  can  affect  the  fre¬ 
quency  distribution  in  the  supercontinuum  generated  in  bulk 
media.  Another  factor  that  strongly  influences  the  intensity 
gradient  is  pulse  self-steepening,  the  effect  of  which  on  the 
supercontinuum  was  discussed  in  Sect.  5.  However,  the  ne¬ 
cessary  condition  for  the  self- steepening  to  be  pronounced 
in  the  bulk  medium  is  self-focusing  of  the  radiation.  Self- 
focusing  is  limited  by  multiphoton  ionization/absorption.  The 
resulting  free-electron  generation  produces  the  blue  shift  de¬ 
scribed  by  the  second  terms  on  the  right-hand  sides  of  (17)  and 
(18).  Thus,  the  plasma  production  accompanies  self-focusing 
in  nearly  all  the  cases  of  high-power  ultra-short  pulse  propa¬ 
gation  in  bulk  media.  Here  we  note  that  in  optical  devices, 
such  as,  for  example,  microstructured  fibers,  the  combination 
of  self- steepening  in  the  Kerr  medium  and  material  dispersion 
of  the  third  and  higher  orders  can  cause  an  asymmetric  blue 
shift  in  the  power  spectrum  without  the  free-electron  gener¬ 
ation  [46].  An  additional  factor  that  is  sometimes  considered 
as  the  possible  reason  for  supercontinuum  generation  is  the 
four-wave-mixing  process  [37, 38].  However,  the  explanation 
based  on  the  four- wave-mixing  process  predicts  that  in  the  red 
side  of  the  spectrum  (co  <  coo)  longer  wavelengths  propagate 
at  larger  angles  to  the  propagation  axis  [37].  The  latter  is  in 
contradiction  with  the  experimental  observations:  no  conical 
emission  was  observed  in  the  infrared  [39].  In  the  numerical 
simulations,  the  results  of  which  are  presented  below,  we  take 
into  account  material  dispersion  up  to  the  second  order  and 
self- steepening  of  the  pulse.  Four- wave  mixing  is  not  involved 
in  our  model. 

The  second  set  of  experiments  in  water  was  performed 
with  the  pulses  generated  by  the  new  laser  system  described 
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in  Sect.  3  (Fig.  la).  The  pulses  were  centered  at  A.o  =  810  nm 
with  a  45-fs  FWHM  duration  and  a  beam  diameter  of  5  mm  at 
1/e2  intensity  level.  The  energy  range  was  from  0.2  to  4  pJ. 
A  series  of  microscope  objectives  was  used  to  vary  the  focal 
length  in  the  experiment.  The  resulting  geometrical  focus¬ 
ing  distances  were  changed  from  10.8  to  73.5  mm.  Numerical 
simulations  in  the  conditions  of  this  experiment  were  per¬ 
formed  according  to  the  system  of  equations  (7)— (1 1)  and  the 
input  pulse  was  in  the  form  of  (14).  The  following  parame¬ 
ters  were  used:  pulse  duration  2to  =  54  fs  (corresponding  to 
45  fs  FWHM),  beam  radius  ao  =  1.7  mm,  geometrical  focus¬ 
ing  distance  Rf  =  16.9  mm,  43.1  mm,  and  73.5  mm,  chirp  pa¬ 
rameter  5  =  0,  the  density  of  neutrals  Aa  =  3.3  x  1022  cm-3, 
band-gap  energy  W  =  6.5  eV,  the  electron  collision  cross  sec¬ 
tion  ac  =  10“ 15  cm2,  and  the  group- velocity  dispersion  co¬ 
efficient  K  =  3x  10  28  s2/cm.  In  the  detailed  experimental 
and  numerical  study  [61],  it  is  shown  that  there  is  competi¬ 
tion  between  laser  pulse  filamentation  and  optical  breakdown 
depending  on  the  geometrical  focusing  condition.  For  a  long 
focal  length,  supercontinuum  generation  can  precede  opti¬ 
cal  breakdown,  while  for  a  short  focal  length  optical  break¬ 
down  precedes  supercontinuum  generation.  There  is  a  critical 
regime  for  the  focal  length  where  both  the  phenomena  can 
coexist  together  in  the  optical  medium.  In  this  regime,  a  ran¬ 
domly  distributed  white-light  beam  is  a  characteristic  sign  for 
the  appearance  of  the  filamentation  and  the  supercontinuum 
generation  near  the  breakdown  plasma  [51,61]. 

In  Fig.  9a  we  present  the  experimentally  obtained  (solid 
line)  and  simulated  (dashed  line)  frequency  spectra  of  the 
pulse  focused  in  water  with  the  geometrical  focal  distance 
of  73.5  mm.  Both  in  the  experiment  and  in  the  simulations 
the  spectrum  was  integrated  over  the  entire  transverse  aper¬ 
ture  and  normalized  to  its  maximum  value.  The  pulse  energy 
was  3  pJ.  Good  qualitative  agreement  between  the  two  spectra 
is  observed.  The  map  of  spectral  wavelengths  corresponding 
to  the  frequency  spectrum  is  shown  in  Fig.  9b.  Similarly  to 
Fig.  4,  the  distribution  of  wavelengths  is  presented  together 
with  the  spatio-temporal  intensity  distribution  /(r,  t)  (equal- 
intensity  contours  are  shown  by  solid  lines).  The  upper  and 
right  insets  demonstrate  intensity  and  phase  distributions  in 
temporal  and  spatial  cross  sections,  respectively.  As  in  air, 
in  water  there  is  also  a  distinct  supercontinuum  source  right 
after  the  self-focused  front  of  the  pulse.  The  location  of  this 
source  is  at  x  =  — 16  fs,  r  =  0.  At  the  same  time,  in  addition 
to  this  first  source,  there  are  multiple  sources  located  in  the 
rings  produced  at  the  trailing  edge  of  the  pulse.  Several  tem¬ 
poral  phase  jumps  shown  in  the  upper  inset  of  Fig.  9b  confirm 
the  formation  of  the  high-frequency  components  in  the  rings. 
The  latter  is  manifested  in  the  conical  emission  rings  observed 
(but  not  studied  in  detail)  in  the  present  experiment  in  water. 
The  radius  of  the  ring  was  increasing  with  the  decrease  in  the 
wavelength. 

The  increase  in  the  amount  of  supercontinuum  sources  in 
water  is  associated  with  1000  times  smaller  critical  power  for 
self-focusing  in  comparison  with  air.  The  location  and  the 
number  of  sources  strongly  depend  on  the  geometry  of  the 
experiment.  By  elongating  the  region  of  filamentation  in  wa¬ 
ter,  one  can  obtain  smooth  growth  of  the  amount  of  the  input 
pulse  energy  redistributed  to  the  white  light.  With  tight  focus¬ 
ing  a  sharp  increase  of  the  white-light  energy  occurs  at  a  very 


FIGURE  9  a  Spectra  obtained  for  the  propagation  of  a  Ti: sapphire  laser 
pulse  with  the  duration  45  fs  FWHM  and  the  energy  of  3  pJ  in  water.  Experi¬ 
ment  {solid  line),  simulations  {dashed  line).  In  the  simulations  the  spectrum 
is  integrated  over  the  whole  transverse  aperture  and  normalized  to  the  max¬ 
imum  value  of  the  spectral  intensity.  Propagation  distance  z  =  —0.66  mm 
from  the  geometrical  focus  position  in  the  water  cell,  b  Spatio-temporal  dis¬ 
tribution  of  the  supercontinuum  sources  corresponding  to  the  simulated  spec¬ 
trum  shown  in  (a).  The  solid  lines  are  the  equal-intensity  contours  plotted  for 
the  following  intensity  values:  In  =  I()4n ,  where  n  =  —  5,  —4,  —  3, . . .  ,  +2, 
/o  ~  1012  W/cm2.  The  upper  inset  corresponds  to  the  temporal  distribution 
of  the  pulse  intensity  {solid  line)  and  phase  {dashed  line)  at  r  =  65  pan; 
the  right  inset  corresponds  to  the  radial  distribution  of  the  intensity  {solid 
line)  and  the  phase  {dashed  line)  at  r  =  0.  Radial  and  temporal  positions 
r  =  65  p,m  and  r  =  0,  respectively,  are  indicated  by  the  white  dashed  lines  on 
the  map 


short  distance;  however,  the  overall  percentage  of  energy  in 
the  spectral  wings  might  be  smaller.  Thus,  the  geometrical 
focusing  condition  of  the  experiment  is  one  of  the  tools  for 
governing  the  supercontinuum  generation.  Another  possible 
tool  is  initial  phase  modulation  of  the  laser  pulse  (see  Sect.  8). 

7  Supercontinuum  from  multiple  filaments 

In  order  to  increase  the  energy  of  the  supercontin¬ 
uum  in  air,  terawatt  pulses  are  used  in  the  experiments.  The 
peak  power  of  such  pulses  exceeds  by  several  hundred  times 
the  critical  power  for  self-focusing  in  air  [18-20].  When  such 
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pulses  are  propagated  in  the  atmospheric  conditions,  a  bunch 
of  filaments  is  formed  in  the  transverse  beam  section.  Multi- 
filamentation  is  a  stochastic  process,  since  the  breakup  of  the 
laser  pulse  into  the  filaments  is  stimulated  by  random  inten¬ 
sity  perturbations  in  the  transverse  beam  section.  The  source 
of  such  perturbations  is  imperfect  initial  beam  quality  [64] 
as  well  as  the  refractive-index  fluctuations  caused  by  atmo¬ 
spheric  turbulence  [65].  Filaments  are  formed  at  different  dis¬ 
tances  from  the  output  of  the  compressor;  their  positions  in  the 
transverse  beam  section  are  random.  Initially,  the  filaments 
are  born  in  the  high-intensity  near-axis  region  and  afterwards 
at  the  periphery  of  the  beam.  With  increasing  distance  from 
the  output  of  the  compressor,  the  probability  to  register  the 
larger  number  of  filaments  increases  [66]. 

Since  the  formation  of  multiple  filaments  results  from  the 
random  intensity  perturbations  in  the  beam  section,  each  fil¬ 
ament  can  be  considered  as  an  independent  source  of  the 
supercontinuum  [67,68].  The  superposition  of  such  sources 
gives  rise  to  a  white-light  laser  pulse  with  comparatively  large 
power.  For  the  propagation  of  a  1 .6-TW,  120-fs  pulse  through 
the  noble  gases,  the  spectral  power  of  1 00- 1 000  MW / nm  was 
registered  in  the  experiment  [13]. 

Study  of  the  supercontinuum  coherence  in  the  course  of 
a  terawatt  pulse  breakup  into  a  bunch  of  random  filaments 
is  important  in  view  of  femtosecond  lidar  applications.  The 
measurements  of  the  coherence  of  the  supercontinuum  pro¬ 
duced  by  the  laser  pulse  focused  into  different  liquids  have 
shown  that  the  relative  length  of  coherence  of  all  blue  spectral 
components  is  the  same  as  the  coherence  length  of  the  main 
laser  wavelength  of  800  nm  [2, 68]. 

In  the  experiments  on  several-filament  formation  from 
laser  pulses  generated  by  a  1-TW,  10-Hz  pulse  of  the 
Ti: sapphire  system  shown  in  Fig.  la,  the  interference  of  ring 
structures  produced  by  two  filaments  was  registered  [64].  The 
pulse  duration  was  45  fs  and  the  energy  14  mJ.  The  regular 
structure  of  concentric  rings  produced  by  the  filament  de¬ 
veloped  earlier  in  the  propagation  is  distorted  in  the  vicinity 
of  the  second  filament,  which  was  formed  later  in  the  propa¬ 
gation  (Fig.  10).  The  maxima  of  the  fluence  in  the  radial 
direction  are  due  to  the  interference  of  the  two  divergent  fields 
of  filaments.  Both  fields  are  at  the  laser  central  frequency  of 
800  nm.  The  interference  pattern  is  the  evidence  for  the  coher¬ 
ence  of  the  two  divergent  fields.  Therefore,  the  sources  of  the 
supercontinuum  blue  wing,  which  are  in  the  rings,  surround¬ 
ing  the  filament,  are  also  coherent. 

In  the  experiments  on  the  observation  of  the  conical  emis¬ 
sion  from  multiple  filaments  we  used  45-fs  pulses,  generated 
by  the  same  laser  system  as  in  the  previous  experiment,  but 
with  a  higher  energy  equal  to  50  mJ.  Fluence  distribution  in 
the  beam  at  the  output  of  the  compressor  was  close  to  Gaus¬ 
sian  with  the  diameter  of  1.9cm  at  1/e  fluence  level.  The 
relative  value  of  the  fluence  perturbations  in  the  beam  profile 
did  not  exceed  5%  and  their  transverse  scale  was  of  the  order 
of  1  mm.  However,  due  to  the  high  peak  power  of  the  pulse 
(Ppeak  =1.1  TW),  the  hot  spots  in  the  beam  profile  were  cre¬ 
ated  already  after  the  1 -cm- thick  CaF2  the  window  located  at 
the  output  of  the  vacuum  pipe  (see  Fig.  la).  These  hot  spots 
developed  into  the  filaments  after  the  propagation  of  4.5  m 
in  the  corridor  air.  With  increasing  propagation  distance  the 
number  of  filaments  as  well  as  the  transverse  region  of  their 


FIGURE  10  Single-shot  transverse  fluence  distribution  in  a  45  fs,  14  mJ 
pulse  at  a  distance  z  =  87  m  from  the  output  window  to  the  corridor  (see 
Fig.  la).  The  number  of  the  neutral-density  filters  in  front  of  the  CCD  cam¬ 
era  is  chosen  so  as  to  saturate  both  the  central  and  the  side  filaments.  The 
interference  of  the  transverse  ring  structures  coming  from  the  two  filaments 
is  observed 

random  formation  was  increasing.  A  detailed  description  of 
this  experiment  will  be  published  elsewhere  [69].  The  follow¬ 
ing  gives  a  brief  account  of  some  essential  results.  Figure  1  la 
shows  the  single-shot  image  taken  by  a  Cohu  4800  CCD 
camera  of  the  transverse  fluence  distribution  at  a  distance 
z  =  30  m  from  the  CaF2  window.  Bright  spots  demonstrate 
the  developed  filaments,  while  pale  spots  show  the  filaments, 
which  are  still  under  development. 

To  register  the  supercontinuum  conical  emission  we  use 
a  Canon  A  40  digital  camera.  The  digital  camera  image  ob¬ 
tained  at  the  same  distance  z  =  30mas  the  CCD  image  is 
shown  in  Fig.  lib.  Colored  conical  emission  rings  are  en¬ 
compassing  the  bunch  of  filaments.  The  radius  of  the  ring 
increases  with  increasing  spectral  frequency  as  in  the  case 
of  a  single  filament.  However,  unlike  the  case  of  the  single¬ 
filament  regime,  there  are  interference  maxima  in  the  radial 
direction.  Note  that  the  interference  pattern  in  Fig.  1  lb  is  reg¬ 
istered  for  the  spectral  components  of  the  supercontinuum  and 
not  for  the  laser  frequency  as  in  Fig.  10.  The  interference  pic¬ 
ture  for  the  bunch  of  filaments  has  a  more  complex  structure 
than  in  the  case  of  the  two-filament  regime.  The  most  contrast 
maxima  are  formed  near  the  most  remote  filament  as  the  re¬ 
sult  of  the  interference  between  the  supercontinuum  formed 
by  this  particular  filament  and  the  rest  of  the  filaments  in  the 
bunch.  With  increasing  propagation  distance  and,  hence,  the 
number  of  filaments,  the  interference  maxima  become  less 
contrasty.  The  interference  pattern  was  changing  from  pulse 
to  pulse  due  to  the  random  process  of  the  filament  formation 
in  the  beam  section. 

On  the  basis  of  these  results  the  following  scenario  of  su¬ 
percontinuum  generation  from  multiple  filaments  can  be  de¬ 
duced.  In  each  pulse  the  supercontinuum  sources  are  created, 
which  are  located  randomly  in  the  transverse  beam  section. 
The  phase  shift  between  the  same  spectral  components  gen¬ 
erated  by  different  sources  depends  on  the  source  position 
in  both  longitudinal  and  transverse  directions.  Since  the  fil¬ 
ament  formation  position  depends  on  the  beam  quality  and 
refractive-index  fluctuations  in  the  atmosphere,  the  phase  shift 
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FIGURE  11  a  Single-shot  transverse  fluence 
distribution  in  a  45 -fs,  50-mJ  pulse  at  a  dis¬ 
tance  z  =  30  m  from  the  output  window  to  the 
corridor.  The  number  of  the  neutral-density  fil¬ 
ters  in  front  of  the  CCD  camera  is  chosen  so 
as  none  of  the  multiple  filaments  is  saturated, 
b  The  digital  camera  image  of  the  transverse 
fluence  distribution  for  the  same  laser  pulse 
parameters  and  propagation  distance  as  in  a 


is  a  random  value  from  shot  to  shot.  However,  during  a  sin¬ 
gle  laser  pulse,  the  supercontinuum  sources  remain  coherent 
since  they  originate  from  the  same  coherent  pulsed  radia¬ 
tion.  Refractive-index  fluctuations  in  the  atmosphere  can  be 
considered  as  frozen  on  the  femtosecond  time  scale.  As  the 
result,  in  the  course  of  multifilamentation  of  each  pulse  we 
observe  the  interference  pattern  from  several  supercontinuum 
sources.  With  propagation  distance  the  interference  pattern 
can  experience  transformations  due  to  the  change  in  the  phase 
shift  between  these  sources.  Besides,  from  pulse  to  pulse  the 
phase  shift  changes  randomly  due  to  the  variations  in  the  field 
of  atmospheric  refractive-index  fluctuations  and  the  intensity 
perturbations  at  the  laser  system  output.  Therefore,  the  inter¬ 
ference  pattern  might  be  different  in  each  single  pulse.  For 
a  train  of  pulses,  the  interference  picture  and,  hence,  the  direc¬ 
tional  diagram  of  the  supercontinuum,  is  averaged  and  has  an 
axially  symmetric  distribution  in  the  form  of  a  narrow  cone. 
The  averaged  power  of  spectral  components  is  proportional  to 
the  number  of  filaments  in  the  beam. 

8  Control  of  the  supercontinuum  generation 

by  means  of  initial  pulse  chirp 

From  the  point  of  view  of  applications,  it  is  import¬ 
ant  to  control  the  nonlinear-optical  transformation  of  a  fem¬ 
tosecond  laser  pulse  into  the  supercontinuum.  In  the  experi¬ 
ment  [19],  performed  in  the  actual  atmospheric  conditions, 
the  supercontinuum  generation  was  controlled  simply  by  ge¬ 
ometric  focusing  of  the  beam.  Initial  wavefront  curvature  to 
a  large  extent  defines  the  filament  starting  position  and  regu¬ 
larizes  the  formation  of  a  bunch  of  filaments  in  the  conditions 
of  turbulence. 

Initially,  the  effect  of  the  pulse  phase  modulation  in  the 
time  domain  (i.e.  temporal  up  and  down  chirps  of  the  pulse) 
on  the  distance  between  the  laser  system  output  and  the  fila¬ 
ment  starting  position  was  experimentally  registered  in  [70]. 
The  results  showed  that,  at  a  constant  pulse  energy,  when  the 
chirp  was  swept  from  a  negative  value  to  a  positive  value,  the 
distance  passed  through  a  minimum  at  the  position  of  zero 
chirp  (Fig.  12).  This  was  explained  in  the  following  way.  With 
either  negative  or  positive  chirp,  the  pulse  has  a  lower  peak 
power.  Hence  it  takes  a  longer  propagation  distance  to  reach 
the  self-focus  or  the  starting  point  of  the  filament.  In  the  case 


of  negative  chirp,  the  compensation  of  the  chirp  by  linear 
dispersion  during  propagation  in  the  neutral  air  also  played 
a  role.  In  the  outdoor  experiments  [71]  the  filament  distance 
was  increased  from  10  to  40  m  by  means  of  the  initial  pulse 
chirp.  Measurements  of  the  supercontinuum  [18]  showed  that 
the  usage  of  the  negatively  chirped  pulses  leads  to  the  increase 
in  the  conversion  efficiency  to  the  near-infrared  part  of  the 
spectrum. 

We  have  performed  a  numerical  study  of  the  filamentation 
and  supercontinuum  generation  control  by  means  of  the  initial 
pulse  chirp.  The  parameters  of  the  radiation  were  X  =  800  nm, 
pulse  energy  =  60  mJ,  duration  of  a  transform-limited  pulse 
at  1  /  e  level  of  intensity  was  ro  =  42  fs,  beam  radius  was 
clq  =  1.5  cm.  The  chirped  pulse  duration  tp(<5)  was  defined 
from  the  condition  of  a  constant  spectral  width  A co  [56].  If  the 
parameter  8  >  0,  the  phase  modulation  is  positive  and  in  the 
medium  with  normal  dispersion,  as  air,  and  the  linear  propa¬ 
gation  regime,  the  duration  of  the  chirped  pulse  increases.  If 
the  parameter  8  <  0,  the  phase  modulation  is  negative,  and  in 
air  the  pulse  duration  decreases  and  then  increases.  Without 
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FIGURE  12  Dependence  of  the  beginning  of  the  filament  on  the  spacing  be¬ 
tween  the  parallel  grating  pair  in  the  compressor,  zm  is  the  distance  from  the 
output  of  the  compressor  to  the  beginning  of  the  filament.  Input  beam  ra¬ 
dius  at  1/  e  intensity  level  ao  =  2. 1  mm;  pulse  energy:  6.2  mJ;  optimum  pulse 
duration:  250  fs.  Positive  chirp  is  at  larger  grating  positions  (towards  the 
right-hand  side  of  the  minimum);  negative  chirp  towards  opposite  positions 
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the  nonlinear  effects,  the  pulse  duration  is  a  minimum  at  a  dis¬ 
tance  equal  to  the  compression  length  Lcomp(<$)  [41]. 

The  conversion  efficiency  of  the  radiation  at  the  laser  fre¬ 
quency  to  the  supercontinuum  is  defined  by  the  whole  process 
of  filamentation  and  depends  on  the  filament  length  as  well  as 
the  presence  of  the  regions  with  high  temporal  gradients  of  the 
electric  field.  The  effect  of  the  initial  pulse  chirp  on  the  fila¬ 
mentation  is  revealed  simultaneously  through  the  two  factors: 
static  and  dynamic. 

The  static  factor  is  in  the  redistribution  of  the  power  in  the 
pulse  at  the  output  of  the  compressor.  Indeed,  for  the  chirped 
pulse,  tp(<$)  >  tp(0)  and  the  peak  power  /peak(<5)  decreases 
with  increasing  absolute  value  of  8  according  to  the  relation: 

/peak(<$)  =  7tall0T0/Tp(S )  .  (19) 


As  the  peak  power  decreases,  the  self-focusing  length 
Lsf(S)  [41]  and,  hence,  the  filament  formation  and  supercon¬ 
tinuum  generation  distance  increase. 

The  dynamic  manifestation  of  the  initial  pulse  chirp  is  in 
the  change  of  the  light-field  localization  in  time  in  the  dis¬ 
persive  medium.  The  character  of  these  changes  depends  on 
the  chirp  sign.  Figure  13  shows  the  fluence,  peak  intensity,  and 
supercontinuum  energy  distributions  along  the  propagation 
axis  z  for  different  values  of  the  chirped  pulse  duration  tp(<5). 
In  each  panel  the  first  plot  corresponds  to  the  fluence  distri¬ 
bution,  while  the  second  plot  shows  the  change  of  the  pulse 
peak  intensity  (thin  solid  line  plotted  against  the  left  vertical 
axis)  and  the  supercontinuum  energy  (thick  solid  line  plot¬ 
ted  against  the  right  vertical  axis).  Note  that  the  fluence  is  the 
value  that  is  usually  measured  in  the  experiments.  One  can  see 
that,  independently  of  the  initial  pulse  chirp,  the  fluence  dis¬ 
tribution  along  the  filament  is  nonmonotonic  and  represents 
the  sequence  of  bright  spots.  In  the  case  of  the  positive  ini¬ 
tial  phase  modulation  ( 8  >  0)  the  spreading  out  of  the  pulse  in 
air  is  stronger  than  in  the  case  of  the  transform-limited  pulse 
(compare  Fig.  13a  and  b).  Simultaneously,  such  localization 
leads  to  the  elongation  of  the  filament  in  comparison  with  the 
transform-limited  pulse. 

In  the  case  of  the  negative  initial  phase  modulation  (8  <  0), 
the  compression  of  the  pulse  takes  place  and  its  effect  on  the 
filamentation  depends  on  the  relation  between  the  compres¬ 
sion  length  Lcomp(<5)  and  the  self-focusing  length  Lsf(<5).  If 
^comp(^)  Lsf(<$)  then  the  pulse  contraction  due  to  the  ma¬ 
terial  dispersion  precedes  self-focusing  and  the  filamentation 
develops  similarly  to  the  case  of  the  transform-limited  pulse. 
If  Lcomp  ( 8 )  ^>>  LSf  (<5)  then  the  effect  of  the  dispersion-provided 
contraction  on  the  filamentation  is  negligible. 

The  optimum  situation  for  the  control  of  the  filamentation 
is  the  comparable  values  of  compression  and  self-focusing 
lengths  with  Lcomp(8)  >  Lsf(<5).  In  this  case  in  the  preliminary 
stretched  pulse  (rp  (<5)  >  200  fs  for  our  parameters),  the  energy 
is  gradually  localized  in  time  in  the  course  of  the  filamen¬ 
tation.  As  the  result,  we  obtain  the  increase  of  the  filament- 
formation  distance  and  the  simultaneous  increase  in  the  length 
of  the  filament  in  comparison  with  the  transform-limited 
pulse  (Fig.  13b  and  c).  Indeed,  in  the  course  of  the  pulse 
temporal  compression,  the  power  in  the  front  slices  becomes 
higher  with  propagation  distance.  A  number  of  slices,  which 
at  the  laser  system  output  did  not  have  enough  power  for  self- 


FIGURE  13  In  a  to  d  the  upper  plot  shows  the  gray-scale  map  of  the  flu¬ 
ence  distribution  in  the  filament  along  the  propagation  axis  z.  The  lower  plot 
shows  the  dependence  of  the  pulse  peak  on-axis  intensity  /peak  ( thin  solid 
line  plotted  against  the  left  vertical  axis )  and  the  supercontinuum  energy  Wsc 
(; thick  solid  line  plotted  against  the  right  vertical  axis )  on  distance  z.  Pulse 
energy  60  mJ,  Iq  ~  1013  W/cm2,  transform  limited  pulse  duration  to  =  21  fs. 
At  z  =  0  the  pulse  has  initial  phase  modulation  a  <5  >  0,  rp(5)  =  100  fs; 
b  8  =  0,  rp(<5)  =  21  fs;  c  8  <  0  (negative  chirp),  rp(5)  =  150  fs;  d  8  <  0  (nega¬ 
tive  chirp),  rp(5)  =  600  fs 


focusing,  accumulate  the  power  P  >  Pcr  at  some  propagation 
distance.  Therefore,  self-focusing  of  these  slices  is  delayed 
in  distance,  and  the  result  of  this  is  the  lengthening  of  the 
filament. 

In  the  simulations  for  the  negatively  chirped  pulse  with 
tp(<$)  ~  150  fs,  the  compression  length  Lcomp(5)  nearly  coin¬ 
cides  with  the  self-focusing  length  Lcomp(<5).  The  localization 
of  energy  in  both  space  and  time  occurs  simultaneously.  In  the 
case  of  such  spatio-temporal  ‘focusing’  the  filament  is  degen¬ 
erated  into  several  bright  maxima  and  its  length  is  decreased 
even  in  comparison  with  the  case  of  the  transform-limited 
pulse  (compare  Fig.  13b  and  c). 

Temporal  compression  of  the  pulse  essentially  affects  the 
spatio-temporal  gradients  of  the  intensity  as  well  as  the  phase 
and,  hence,  the  conversion  efficiency  of  the  initial  laser  pulse 
to  the  supercontinuum.  This  is  confirmed  by  the  simulation 
results  for  the  peak  intensity 


W (r  =  0,z)  =  max  {I(r  =  0,  z,  r)}  ,  (20) 
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and  the  energy  contained  on  the  blue  side  of  the  supercontin¬ 
uum 

■^•max 

Wsc(z)  =  j  S(X,z) dX,  (21) 

^min 

where  Xm[n  =  500  nm  and  A.max  =  700  nm.  As  mentioned  ear¬ 
lier,  the  values  (20)  and  (21)  are  shown  in  the  second  plots  of 
panels  (a)-(d)  in  Fig.  13. 

By  following  the  dependence  of  the  supercontinuum  en¬ 
ergy  on  distance,  one  can  notice  that  for  all  four  values  of  the 
chirped  pulse  duration  rp  (5)  the  growth  of  the  supercontinuum 
energy  Wsc  occurs  at  the  same  values  of  z  as  the  jumps  in  the 
peak  on-axis  intensity  /pea k- 

In  the  transform-limited  pulse  the  change  of  the  peak  in¬ 
tensity  along  the  filament  is  comparatively  smooth.  The  lead¬ 
ing  fronts  of  the  intensity  are  not  very  steep  and  the  corres¬ 
ponding  phase  gradients  are  moderate.  The  resulting  conver¬ 
sion  efficiency  to  the  supercontinuum  is  low  and  the  energy 
contained  on  the  blue  side  of  the  supercontinuum  Wsc  does  not 
exceed  0.01%  of  the  initial  pulse  energy  Wo  (Fig.  13b). 

In  the  case  of  the  positive  phase  modulation  the  intensity 
peaks  at  the  end  of  the  filament  lead  to  the  increase  of  the  su¬ 
percontinuum  energy  up  to  1%  of  the  initial  pulse  energy  Wo 
(Fig.  13a). 

In  the  case  of  the  negative  phase  modulation  the  filament 
consists  of  a  number  of  contrast  intensity  maxima.  The  cor¬ 
responding  temporal  gradients  of  the  intensity  and  the  phase 
are  large  enough  to  lead  to  the  increase  in  the  supercontinuum 
energy  up  to  3%  of  Wo  (Fig.  13c  and  d).  The  highest  super¬ 
continuum  energy  is  achieved  in  the  shortest  filament  with 
the  chirped  pulse  duration  rp(<5)  =  150  fs  (Fig.  13c).  This  fil¬ 
ament  consists  of  the  three  separated  regions  with  high  peak 
intensity.  In  these  regions  high  localization  of  the  radiation 
provided  due  to  both  temporal  compression  and  self-focusing 
is  replaced  by  the  plasma-induced  defocusing.  The  resulting 
phase  gradients  lead  to  the  effective  conversion  of  the  laser 
pulse  energy  into  the  supercontinuum  energy. 

Thus,  in  the  case  of  the  filamentation  of  the  negatively 
chirped  pulses,  the  efficiency  in  the  generation  of  the  high- 
frequency  components  of  the  supercontinuum  is  much  higher 
than  in  the  case  of  the  transform-limited  pulse.  The  super¬ 
continuum  energy  increases  by  two  orders  of  magnitude  in 
a  wide  range  of  the  negatively  chirped  pulse  duration  rp(<5) 
from  tp(<5)  =  100  fs  to  tp(<5)  =  800  fs. 

Note  that  the  energy  contained  on  the  red  side  of  the  su¬ 
percontinuum  is  weakly  dependent  on  the  initial  phase  mod¬ 
ulation  of  the  laser  pulse.  This  is  due  to  the  fact  that  the  long- 
wavelength  components  of  the  supercontinuum  are  generated 
early  in  the  front  of  the  pulse.  The  steepness  of  this  particular 
region  is  nearly  not  changing  in  the  case  of  the  filamentation 
of  initially  chirped  pulses. 

9  Conclusions 

We  have  performed  an  experimental  and  numeri¬ 
cal  study  of  the  nonlinear-optical  transformation  of  a  fem¬ 
tosecond  laser  pulse  in  condensed  and  gaseous  media  by 
the  example  of  water  and  atmospheric  air,  respectively.  We 


have  shown  that  the  spatio-temporal  modification  of  the  pulse 
shape  and  the  corresponding  frequency  angular  modifica¬ 
tion  of  the  pulse  spectrum  result  from  the  same  nonlinear- 
interaction  of  the  laser  radiation  with  the  medium.  This  in¬ 
teraction  is  due  to  the  joint  effect  of  the  Kerr  nonlinearity 
of  neutrals  and  the  nonlinearity  of  the  laser-produced  plasma 
as  well  as  such  propagation  effects  as  dispersion,  diffraction, 
and  self-steepening  of  the  pulse  fronts.  Formation  of  filaments 
with  high  spatio-temporal  localization  of  energy  provides  us 
with  the  large  nonlinear  interaction  length,  which  leads  to 
a  dramatic  change  of  the  pulse  shape  and  transformation  of  the 
pulse  spectrum  to  the  supercontinuum. 

Transformation  of  a  narrow  laser  pulse  spectrum  into  the 
continuum  is  due  to  the  self-phase  modulation  of  the  light 
field  in  both  space  and  time.  The  temporal  phase  gradient  is 
responsible  for  the  new  frequency  components,  while  the  spa¬ 
tial  gradient  is  responsible  for  the  angular  divergence  of  these 
newborn  frequencies.  As  a  result,  the  supercontinuum  cone 
emission  is  observed.  High  spatio-temporal  localization  of  the 
radiation,  which  is  observed  in  the  experiments,  results  from 
the  dynamic  process  of  the  replenishment  of  energy  in  the 
filament.  The  replenishment  process  encompasses  a  signifi¬ 
cant  transverse  area,  which  is  larger  than  the  beam  diameter 
at  the  laser  system  output.  In  this  region  the  dynamic  ring 
structure  is  formed  in  the  intensity  distribution.  The  correct 
numerical  consideration  of  the  low-intensity  background  at 
the  beam  periphery  allowed  us  to  obtain  satisfactory  quanti¬ 
tative  agreement  between  the  simulated  and  experimentally 
registered  supercontinuum  conical  emission  angles.  Both  the 
experiment  and  the  simulations  were  performed  for  the  propa¬ 
gation  of  a  250-fs,  10-mJ  pulse  with  a  central  wavelength  of 
800  nm  and  a  beam  radius  of  3.5  mm  at  1/e  intensity  level. 

We  have  shown  that  the  sources  of  the  supercontinuum 
blue  wing  are  in  the  ring  structure  surrounding  the  fila¬ 
ment.  Another  contribution  to  the  high-frequency  spectral 
components  comes  from  the  trailing  edge  of  the  pulse,  where 
the  positive  frequency  shift  originates  due  to  the  shock  forma¬ 
tion  at  the  back  of  the  pulse.  This  source  of  the  supercontin¬ 
uum  does  not  reveal  the  increase  of  the  divergence  angle  with 
increasing  spectral  frequency.  The  newborn  high  frequencies 
mainly  propagate  along  the  beam  axis. 

In  water  we  experimentally  obtained  two  characteristic 
features  of  the  spectral  broadening.  At  low  laser  energies 
(0.4- 1.8  |xJ)  the  spectra  of  a  170-fs  pulse  show  a  symmetri¬ 
cal  broadening  towards  lower  and  higher  frequencies,  while 
at  higher  pulse  energies  (2-200  |xJ)  the  typical  long  pedestal 
emerges  on  the  blue  side.  The  minimum  wavelength  of  the 
blue  spectral  wing  (~  400  nm)  remains  almost  constant  for 
an  increase  in  the  pulse  energy  by  two  orders  of  magnitude. 
This  is  the  evidence  for  the  peak  intensity  clamping  inside 
the  filaments  produced  by  an  ultra-short  pulse  in  water.  Nu¬ 
merical  simulations  of  a  45-fs,  3-|xJ  pulse  propagation  in 
water  showed  good  qualitative  agreement  between  the  simu¬ 
lated  and  experimentally  obtained  frequency  spectra.  Multi¬ 
ple  sources  of  the  high-frequency  components  of  the  super¬ 
continuum  are  found  to  be  located  in  the  rings  formed  due  to 
the  defocusing  of  the  pulse  in  the  plasma,  produced  due  to 
both  multiphoton  ionization  and  the  avalanche. 

In  the  experiments  on  multiple  filamentation  of  a  45-fs 
pulse  with  the  energy  of  14-50  mJ  and  the  central  wavelength 
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of  810  nm  we  have  registered  the  phenomenon  of  the  inter¬ 
ference  of  the  conical  emission  spectral  components,  gener¬ 
ated  by  different  filaments  from  the  same  bunch  of  filaments. 
Therefore,  the  white-light  laser  pulses  produced  by  multi¬ 
ple  filaments  are  coherent  in  spite  of  the  random  character 
of  multifilamentation  associated  with  the  beam  quality  and 
refractive-index  fluctuations  in  the  turbulent  atmosphere. 

The  conversion  efficiency  of  the  radiation  to  the  supercon¬ 
tinuum  is  defined  by  the  whole  process  of  filamentation  and 
depends  on  the  filament  length  as  well  as  the  presence  of  the 
regions  with  high  temporal  gradients  of  the  light  field.  If  the 
initial  pulse  is  negatively  chirped,  the  efficiency  in  the  gener¬ 
ation  of  the  high-frequency  components  of  the  supercontin¬ 
uum  is  much  higher  than  in  the  case  of  the  transform-limited 
pulse.  Numerical  simulations  of  the  propagation  of  a  42-fs, 
60-mJ  pulse  in  air  showed  that  the  supercontinuum  energy  can 
be  increased  by  more  than  an  order  of  magnitude  in  the  case  of 
the  negatively  chirped  pulses. 
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